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PREFACE 


Tne author has attempted to illustrate in a somewhat 
brief manner the application of certain statistical techniques 
to the analysis of core sampling data. The statistical areas 
of frequency distributions, analysis of variances, and to a 
lesser degree, sampling, provide the basis for the study. 

A convenient reference system is used in the thesis. All 
equations are numbered as (Chapter, Number). Thus, equation 
(2-4) signifies the fourth equation of the second chapter. 

The notes are listed together before the bibliography and 

are numbered consecutively within each chapter. All mathemati- 
cal notations and significant terms used are listed and 

defined in Appendix A for ready — 

Actual field data was obtained for the study. The 
author is indebted to the following people who provided the 
fleld data upon which the project was based: Mr. Mack C. 

Colt and Мг. Wendell Weatherby of Iola, Kansas; Mr. Schermerhorn 
of Tulsa, Oklahoma; Mr. Ray Plummer of Chanute, Kansas; and 

Mr. Carl Pate and the Oil Field Research Laboratory of 

Chanute, Kansas. 

The author wishes to express his appreciation to Dr. 
Charles F. Weinaug for his overall direction of the graduate 
program which led to the thesis, and to the Bureau of Supplies 
and Accounts, United States Navy, whose sponsorship made the 
thesis possible. 

The author is most deeply indebted to Dr. Floyd Preston 
whose close personal guidance and wise counsel were instru- 


mental in bringing the thesis to completion. 
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CHAPTER I 
INTRODUCTION 


Background 

Statistics may be considered in two senses. One con- 
cepulon of statistics is that of a collection of numerical 
or quantitative data, i.e., figure data, such as numerical 
data on births or unemployment. Statistics in the second 
sense is less well known, and refers to the techniques of 
analyzing data for decision-making. It may be thought of as 
the science of decision-making in the face of uncertainty. 

The application of statistics in this second sense to 
petroleum engineering problems is the purpose of this thesis. 
The petroleum engineer by the very nature of the realm in 
which he must operate is continually faced with making de- 
cisions based upon fragmentary and inconclusive information 
concerning the sub-surface of the earth. Because of the 
extreme complexity of the geological processes of erosion, 
material transport, deposition and burial of material, the 
porous media forming petroleum reservoirs are extremely 
heterogeneous. Physical properties can vary extensively from 
place to place within individual reservoirs. The extreme 
cost of sampling these reservoirs forces the engineer to 
construct mental and mathematical models from fragmentary 
information. Mathematical statistics would seem to provide 
a valuable Cool {or such model construction. The utilization 


of the science of statistics to provide systematic analysis 





techniques to data that is available may provide additional 
relevant information as an assist to the decision-making process. 

It is the purpose of this thesis to test the applicability 
of certain statistical techniques for creating mathematical 
models of property variation within petroleum reservoirs and 
to show how certain statistical techniques can be used to 
extend the interpretation of fragmentary data. The study is 
limited to the analysis of core samples. Statistical techni- 
ques are applied to actual field data obtained from five 
fields. Specific results are presented to demonstrate calcula- 
tional techniques. 

me cirst approach in the study is to examine the 
probable frequency distributions of the various properties, 
testing to see if property distributions within individual 
fields satisfy the Gaussian normal distribution and if not, 
to describe the distributions by a family of generalized 
frequency curves known as the Pearson system of frequency 
curves. An initial study of the distributions is considered 
essential, prior to the application of the other statistical 
Techniques such as sampling and analysis of variances, since 
these depend upon the nature of the frequency distribution 


of the data under study. 


Origin of Statistics = 


The origin of the modern science of statistics may be 
traced to the mid-seventeenth century when two astute French 


mathematicians, Pascal and Fermot, were presented with a 





problem involving a game of chance and the interpretation 
of the probabilities associated with it. Their work led to 
solutions, not only of the problems proposed, but of more 
general ones. The methods employed by Paseal may be said 
to represent the beginning of the mathematics of probability, 
about which modern statistical theory centers today. The 
publication by Laplace in 1812 of the epoch-making "Theorie 
Analytigue des Probabilities" laid a firm foundation for 
this theory.” 
Statisticians in general from Pascal onward sought a 
method of describing the nature of the distributions of 
pucucecerbeobs. Table pS shows the most notable of the 
methods developed to portray chance effects from the time 
of Pascal until Karl Pearson set forth his system of generalized 
frequency curves in 1895. 
The usual purpose of frequency distributions is to 
represent a sample of actual data drawn from a much larger 
or even infinite population. Even though a sample may be 
composed of a relatively small finite number of observa- 
tions, it may be reasonably representative of the larger 
universe from which it was drawn. Since it is virtually 
impossible to measure all the items comprising a universe, 
it is necessary to form a notion of the larger group from 
the study of a sample. Thus by constructing a hypothetical 
infinite population of which the actual data is regarded as 
constituting a random sample, an understanding of the law 


O distribution of chance effects of the hypothetical 
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population may be obtained and specified by a few parameters. 

For example, to gain an understanding of the porosity 
characteristics of an area, it may be considered that the 
earth contains an infinite number SS measurements 
where it is impossible to measure each and every one of the 
porosities. By fitting a curve to a frequency distribution 
of the actual data obtained from a core sample, it is 
attempted to describe what appears to be the general form 
of the curve for the entire porosity population. 

For statistical work the normal or Gaussian curve is 
probably the best known and most heavily relied upon 
frequency distribution of those shown in Table I, particularly 
in the theory of sampling. However, this distribution 
function does not apply to skewed frequency distributions, 
although numerous sets of data defy the normal curve and 
exhibit markedly skewed distributions. The Gaussian school 
of statisticians regarded skewness as a by-product of 
sampling and believed that skewness could be made to disappear 
completely if an infinite number of observations were available. 

With the recognition that the normal curve was not 
sufficient to characterize all natural observations, it 
became apparent that it was necessary either to devise 
methods of describing the most conspicuous departures from 
the normal distribution or to devise generalized frequency 
curves to describe distributions as they actually exist in 
the observational sphere. 


Karl Pearson followed the latter course and showed that 





a set of frequency curves could be obtained by assigning 
values to the parameters in a certain first order differential 
equation which has its basis in the theory of probability. 
This approach is covered in Chapter II and the application 


of these curves to the actual field data is shown in Chapter III. 


Core Samples 


The taking of core samples from a reservoir has been an 
accepted practice for the past hundred years. At first, well 
Samples and cores:had but one purpose--to locate 011.2 Thus 
it was necessary to take a sample from every well drilled. 
Today, however, with the advent of other techniques such as 
electric logging to aid in locating the oil the primary 
reason for taking samples has shifted to serve as a source 
of information about the reservoir and its contents. At 
present, core samples provide certain numerical parameters 
by which the field may be described. The most common is 
the arithmetic or weighted average of the various properties. 
Ihe range and variances of the distribution of the properties 
are other easily obtained and useful parameters describing 
a field. 

A simple histogram showing the distribution presents 
visually the characteristics of a field. “А cumulative 
frequency curve on graph paper will present the distribution 
and allow easy determination of such parameters as the median 
and possibly the mode. From a cumulative frequency curve an 


estimate of the percentage of a distribution which is above 





a specified minimum point may be obtained. The range of a 
variable within any set quartile is likewise easily deter- 
mined. In summary, graphical methods of statistical data 
presentation permit certain numerical parameters to be 
obtained with relative ease. 

Since it is no longer necessary to take samples from 
Every well> it is desired to determine the number of wells 
that should be core-sampled to provide information needed 
E itn an acceptable probability of obtaining reliable results. 
At present the number of core samples to be taken is deter- 
mined somewhat intuitively with a wide variation of opinion 
as to what is the necessary number. Are there statistical 
techniques available to serve as a guideline in determining 
how much information is needed and how such data should be 
interpreted? 

The method of interpretation of core data may be paramount 
Since raw data by itself, irregardless of the amount available, 
often supplies much information that is irrelevant and 
immaterial. It is the object of statistical processes 
employed to exclude this irrelevant information and to 
isolate the whole of the relevant information contained in 


ри data. 


Data and Techniques 


The science of statistics consists of (1) collecting, 
(2) presenting, and (3) analyzing quantitative data. The 
data used for this study consist of the physical properties 





of oil fields, namely permeability (K), porosity ($), 911 
saturation (Sy), and water saturation (S )> as obtained from 
core samples Taken from five different fields. Table 17" 
shows the type and amount of data considered. The five fields 
will be referred to as Field 1, Field 2, Field 3, Field 4, 

and Field 5. The numbering system has no significance other 
than the fields being numbered in sequence as data were 
obtained for this study. 

Pields 1, 2, 4, and 5 are located in southeast Kansas 
le Hield 5 is in northeast Oklahoma. Field 1 IS in 
County "A" while Fields 2, 4, and 5 are in County "B'. 

The core samples for all five fields were taken and 
analyzed by the same laboratory with the same coring and 
analyzing techniques used for all fields. Thus even though 
there may have been errors made in arriving at the absolute 
values of measurements of the different properties, especially 
With regard to the fluid saturations, the errors may be 
considered consistent, and a relative comparison of the 


data may be made with a certain degree of confidence. 


TABLE II Summary of Core Sampling Data 





No. of Wells No. of X No. of @ No. of S No mor 


Field ү] 





samples Samples samples samples 
1 LO 1503 794 794 794 
2 14 630 345 345 205 
> 7 195 l 152 152 
Y 19 = 316 316 201 
5 AOL - 162 1672 1129 





Standard coring techniques were employed in taking the 
vertical cores of approximately 20 feet in length.  Measure- 
ments of each property were made approximately every Six 
inches in the oil productive section with approximately 
twenty samples being obtained for each well. The depth to 
the pay zone for the different fields varied somewhat, but 
they all could be considered as shallow fields with pay 
zones at depths between six hundred and a thousand feet. 

The well spacing in most instances was approximately 
four hundred feet. In each field more than 50 percent of 
the wells were cored, with coring data available for all 
properties with the exception of permeabilities for Fields 
4 and 5. | 

A map of each of the first three fields showing the 
locations of the wells cored is given in Appendix C. 

The above mentioned data were used for the following 
statistical studies which constitute both the method of and 
scales Lon Гог this thesis: 

l. Analytical Fitting of Data to the Pearson Generalized 
Frequency Curves. The study includes the technique for 
selecting the appropriate Pearson type curve, for 
fitting the data to the selected curve and for measuring 
the goodness of fit of the data to the curve. In addition 
the data were fitted to the normal or Gaussian Curve 
and its goodness of fit determined.  Permeability data 
were converted to logrithms and the resultant distribu- 


tions were analyzed by the above methods. 





10 
2. An Analysis of Variance Study of Well Property Means 
and Variances. 
5. Use of Certain Sampling Methods as a Way of Estimating 
Certain Population Parameters from Point Estimates. 
Пе Application of Additional Statistical Techniques To 


Core Analysis Are Briefly Discussed. 


Statistics, in its many ramifications, is an exceedingly 
complex subject and much too involved to be thoroughly 
covered in a paper of this nature. The techniques presented 
are by no means the only methods of analysis available. In 
order to permit the reader without a thorough background in 
theoretical statistics to gain an understanding of the 
material presented, the discussion of the theoretical proofs 


and principles involved have been held to a minimum. 





CHAPTER TI 
FREQUENCY CURVES 


Introduction 

One of the most important practical problems in mathe- 
matical statistics is the obtaining of a relatively simple 
yet accurate representation of the frequency distribution of 
any set of data under consideration. Some knowledge of the 
frequency distribution of a set of data should be obtained 
before a statistical analysis is attempted. 

There are essentially three methods of describing 
frequency distributions of one variable; namely, the grephi- 
cal method, the method of averages and dispersions, and the 
method of theoretical frequency functions or curves. These 
three methods of describing frequency curves will be briefiy 
compared to their pelative merits. 

The graphical method allows a large amount of data to 
be condensed to an easily presentable form. An inherent 
weakness of this method is the inability to quantitatively 
compare distributions of two or more sets of data. One may 
State that two distributions are somewhat the same or that 
they are somewhat different, but the degree of sameness or 
difference cannot be quantified by observation alone. This 
lack of numerical description of the distribution by the 
graphical method precludes its use as a comparing method 
except in the most elementary studies. 

Figure 1 shows two histograms, one representing the 


frequency distribution of the porosity of Field 1 and the 
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The second method, involving the use of averages and 
dispersion, does give a numerical description of the data 
in Figure 1, but it does not give a functional relation 
between the values of the variable X and the corresponding 
frequencies. The numerical description in terms of X and 


o where: 


| Field X б 
X = Arithmetic mean 1 19. 4% ШОО; 
с = Standard deviation 2 19.7% qr 


This numerical description would further indicate that the 
distributions are nearly the same which again would be some- 
what misleading. 

The third method, an analytical method of describing 
irequency distributions shows that Field 1 is a generalized 
frequency distribution of the Pearson Type lg: whereas Field 
2 is of the Pearson Type iV, and may be described by the 
"d = .725, ô = -.294, and а, 


pespectively. This method indicates that the distributions 


parameters a = 1.86, 8 2.11, 
of the two fields are not similar which is contrary to what 
one would assume by use of the first two methods of describ- 
ing frequency distributions. 

The reader should not be unduly concerned at this point 
as to how the values of а, and 6 were determined and what 


Type ig and IV, mean as the following sections will present 
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ENFEhcPOuEMnMdrPSeussron of Pearson's generalized frequency Curves. 


Pearson's Generalized Freguency Curves 


After it was recognized that the Gaussian or normal 
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eunve tallea Lo describe the distribution of many of the 
observed data, Pearson proceeded to develop generalized 
frequency curves that could characterize the various types 
of unimodal frequency distributions encountered. 

In deciding on a system of curves for describing 


frequency distributions, it was realized that: 


1. Any expression used should be a graduation formula, i.e., 
it must remove the roughness of the data. 
2. An expression, should not involve too many high moments 
to calculate constants, for thereby accuracy is reduced. 
5. There should be a systematic method of analysis applicable 


to all possible types of frequency distributions. 


Then, considering the most obvious characteristics of frequency 

distributions, it may be considered that they generally start 

at zero, rise to maximum, and then fall at the same or often 

at a different rate. At the end of the distribution there is 

often high contact. Mathematically, the above implies that 

a series of equations Y = f(x) or Y = f(t), must be chosen 

so that in each equation of the series dy/dx or dy/dt = 0 for 

certain values of x or t, namely at the maximum and at the 

end of the curve where there is contact with the axis of x or t. 
The above suggests that the frequency function may be 


represented as a solution of the differential equation: 
Ё _ ду _ У(а-% 
(2-1) at 


Since: 


a. For a value of t, t = a, dy/dt = 0 and the required 
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maximum is obtained and, 
b. As Y approaches zero, the derivative dy/dt also approaches 


zero thus giving contact at one end of the curve. 


Assuming that f(t) may be expanded in a converging 


power series, equation (2-1) may be written: 
IRA a-t 

(2-2) EE: Е NE E IL СЫ 
Xy di b. + 19 b, +° + vere 


where the mean of the distribution is taken as the origin, 
and the abscissae are measured in units of the standard 


deviation such that: 


(2-3) E. discancellrom origi 


O 


Tne above suggests that significant frequency functions 


Y = f(t) may be found among the solutions of (2-2) subject to 


the following restrictions: © 
EN 
(2-4) و“‎ = f(t) dat = 
b 
2 
b4 | 
(2-5) a, =f t f(t) dt 20 
b 
2 
Dear 
(2-6) an = f F(t, dt 
Do 
where: 
Han uli 
(2-T) о ==] t Ооа 
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and, 
N М. 
(2-Ta) os E C = E 
Vu 2 Ho 
2 
See Chapter III for definition of the moments Ho» uo > e... Ho. 


Clearing (2-2) of fractions, multiplying through by t апа 
integrating over the range r to s (where r and s are the 


extremes of the range of variation for t) with respect to t 


5 


gives: 


S 
(2-8) [aft"vat-», ftPay-p, Jt ay-p, [s Pay- [o үа ] = о 
eG 


But through integration by parts: 


S S 
(2-8a) ty] - (t? x - nfe?^ v at ] 
Der t-p 


and if the frequency funetion, when multiplied by t? vanishes 


at the limits of the distribution, r and s, we have s 


n 
оо) f t dy = -n N а 


2 


which leads to the recursion formula for moments: 


(2-9a) a, ain a, by + (1+1) ab, + (n-2) a b, =a 


n+] 72 n+1 


Giving n successively the values O, 1, 2 ..., from (2-8) and 


(2-9) noting that N the total frequency cancels out and 


letting Qo = meg = l, one EE 


1 = О, Q 


2 
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а + b. E. TO 
bo + © o EIS 
(2-10) а + 5 bi + 4 az Da *o = 0 
а а + 206 7 79501 + 5 a b, + .... = 0) 
a ЕЕ petu ab, + (n+2) ar, Bo + = Q] 


Assuming that f(t) converges so rapidly that terms involving 
the third and higher powers of t may be neglected, a simul- 


taneous solution of (2-10) yields: 


= D = 23 
2t1+20)” 2. 2(1+26) 


(2-11) 
енко b. = Ô 
о 2(1+26)? 2  2(1-26) 
where, / 
аре Е о а, - 6 
The value of a, where: 
a 


E a = $0528] 


represents the distance between the mean and the mode, which 
is defined by Pearson as the skewness of the distribution. 


Thus from the differential equation: 


dy Na 12 a-t 


(2-2) = > 
Y dt (с by + by t+b,t 


INE 





which has its basis in the theory of probability, the para- 


b.,: and b, have been determined in terms of 


O |. 2 
the moments and expressed in terms of © а Па ES 


meters a, b 
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transformation is desirable to permit the use of the (а, 6) 
chart for identification of the appropriate type curve that | 
fits a set of data. With the above parameters defined, 
Pearson's family of generalized frequency curves may be 
developed. 

In the family of curves there are three main types, 
nine transitional types which are special conditions of the 


three main types and the normal curve. 


Integration of the Differential Equation 


For 6 £0, b, 40, the denominator b, ^ b, t + b, t^ is 


a quadratic which can be written in the form b (t - гу) (© - Po). 


Mus, 
(2-11) i dy _ اھ‎ аи a MU -- е 

E + b; DD b, (t - г1)(% - Po) 
where: 

- D- + /b 2 -1Ъ Ъ - b. - VD - 4b.b 

(2-15) r_ = NES D ee ron and T = ل‎ 1 29 

т 2 р 2 2 р 

2 2 

Upon substitution of az and ó in (2-15), for a, bo» b4, 8nd b: 


У 15(5 = Б 
(2-15a) D. ng зм ымы „ar 


where: 


2” - 146(6-2) and 


o eT a E 
o > 5 (6+2) В 9 JD 


e 2 2 5 25 


р = а 
(2-15b) 
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Also, by the method of partial fractions," 











A 
(2-16) RET O = E ع‎ + B ] where 
b,(t-r, ) (t-r,) b, t-r,  t-r, 
а-г a-r 
@-17) p o. and ae = 
Tes pore i 


Hence from (2-14): 
I ° | 
(2-18) y dy Bu t-r4 г> DoD Kr) 


Upon integration of (2-18): 











ah а-г- 1 а-г, 
(2-19) Log Y = bs (Flor, =.) Log(t-r, ) + Dp orp log(t-r,)+log C 
Hence: 1 (851 ; Lo 
ba ‘ra -r Do ы =]? 
(2-20) ү = с(є-р.) 2 1 2 (t-r,) 2 2 1 
1 2 

( al (2772 al ma 
2-21 Let m. = == and m. = — 

l b. га 2 2 b. Put] 
Thus (2-20) reduces to: 

12 


o(t-r,) (t-r,) 


(2-22) » 27 


Upon substitution of а апа ó in (2-21) 


ô WOS 
CE 


1426 
5 E and 


il 


ا 
| )2-23( 


|| 


m 


> = - ees. E 


For -4 < ô < O, the r's are real and opposite in sign; for 


що, ава Б < 46(6+2), the r's are complex; and for ó > O 
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and 45(6+2) < Os the r's are real and of the Same sign. 
these three conditions with the additional condition that 
о, A O establish the criteria for the three main types of 
Pearson's frequency functions designated Type I, IV, and 
um. The boundaries of these areas, the curve (2+36) u = 
у (1426) (2-6) which intersects the Type I and Type VI areas 
and the line 6 = - 1/2 contains the points which correspond 
to the transitional types. The numbering system of the main 
les, i.e., I, IV, VI, is that established by Pearson and 
ls used in tnis study to provide a standard reference to 


other literature on the family of curves. 


Analysis of Data 


For the purposes of analyzing the field data under 
consideration, it will be shown that Pearson's three main 
types of curves: Type I, Type IV, Type VI, plus the 
transitional Type III and the normal curve will suffice to 
describe the data. Therefore, the development of the other 
transitional types will not be described in detail. Only 


10 


the equations and conditions are shown in Table TER On 


the following pages derivations are given Tor the five types 


of frequency distribution functions used in This thesis. 


туре т 
When the E'S in (2-14) are opposite in sign, (2-22) is 


written as: 


(2-21) ج‎ o(t-r,) } (2-6) 2 
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Equation (2-24) is called Type I of the Pearson system. The 


conditions on a and Ó are: 
од # 0, -1 < 6 < 0 (6 А - 1/2), (2135) a,” £ 4(1+26)° (2+6) 


C is determined by setting the area i.e., the total probability 


equal to unity or, 





“2 ii Па 
(2-25) Bee) l: i) dt = 1 
r a 
1 
Substituting: 
6) t-r- 
2-2 W = 
Pur, 
ETT ME MEE m 
(2-27) Ern S (1-0) ^ aw - 1 
О 
Hence from the definition of the Beta function: +* 
ту + m, +1 
(2-28) Sn) B (m. de o J); = l 


Thus C may be expressed either in terms of the Beta function, 


or alternately in terms of the Gamma function, through use 


of an identity given by Whitaker and Watson: ^2 


| | м 
со) C = m +m,+l š 


(rom) ^ ©  B(my+l,Myr1) 


IF (m, +m,+2) 


тут, 
(22-2) r(m+1)T(m,+1) 


The other parameters are: 
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a, + SD ت‎ 
. 


ЕЕ зз Us 


146 E^ 
ur 2) 


, у= = 


ш Еј = Че р -1 
fora, > 0, гі <0< г, ала EN < 1221. 


The range of the curve is (mj; Po). The curve will be U- 
shaped if both m's are < O, J-shaped if the m's are opposite 


in sign, and bell-shaped if both m's are > o, 14 


Type IV 


Conditions: a ҒО” 0 5 0. = < 4(6+2). The condi- 
tions imply that the гізіп (2-14) are complex and therefore 
the second main Type of curve, Type IV, may be determined. 


Thus (2-15) ean be written: 


T4 = 5 =-r + iS where r = 5 and 8S = 5 


-% -i Y-D 
Pen ЫШ. 
E Ic VI 
EN r2 ып loco с ШІ 
where: 
a 
~2 (352) 2 а 1--26 
/ -D 
THUS: 


| m 
(2-22) Y -C(t-m) ^ (t-r becomes 


2) 
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va 
| -m 
5 E. e 2 t+r-iS 
(2-30) Yt? = 1 (5+2) + 32] (ERES) 
and lo 
сі 
a-bi, = ¿e tan”? b E KC - tan ™ а/о) 
a+bi ü ji 
the frequency function can be written: 
-m -1 /t+r we 
(2-533) v=Cl(ter)?+s2) — QV tan (7) . 2 


To determine C on setting the area of the curve over the 


interval (-o, e), equal to unity, Craig"! shows that: 
2m-1 
% Some 5 
(2-52) ME UICE 
where: 
(2-33 ) G(2m -2, v ) =f "san (2m -2) со Рао 


О 
The term d is defined as: 
е-е 
(2-33a) 6 - (Ë - tan" 5) 


The function G(2m -2, v ) is obtainable in tabular form from 
Pearson, 19 | 

For this thesis, the function G(2m -2, v ) was generated 
by à special numerical integration procedure involving 


Í; 


Gaussian coefficients. 


Type VI 


Conditions: az оо 2, > 46(6+2). The 


conditions imply that the r's in (2-14) are real and of the 
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same sign, thus the third main type of curvo, Туре ү с 


obtained. 
Tg TA 
(2-22) Y = C(t-r-) (t-ran) 
1 2 
An alternative simplified form is: 
m m 
(2-31) у = 0(Z °)(Z - a) Í 


Craig! Shows that: 


T2597 + Bien 
(2-35) C one аси 


P(m +1) Г(-ш, -m4 -1) а 


where: 
(2-36) 2 = t-r, 
(2-37) Q = Py PS 


Range of curve is (ry; c). 


Normal Type Curve 


Conditions: а = ô =O, The original differential 
equation: 
Cc.) і ау _ en 5 
y q Bo + b. t + 55 t 
reduces to: 
lay... 
when: 
a, = 6 = 0, since from equation (2-11) 





== 





2 
Непсе: 
(2-58) Y z -t dt 
ү 
Upon integration: 
+2 
(2-29) ОБ У = = = + 105 С 
ES 
(2-40) | Усе? 
where: 
¿20 __М 
џет о 


where N = total number of observations. 


Type III 


Conditions: Az ХО, 6 =0. From equation (2-11) for 


E > b. = 72 
2 
bees d b, = 0 


Y dt b tb рт 
becomes: 
Q 
жиш 
(2-41) M e ee 
y dt | 


which yields after integration: + 
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— ل‎ -2+ 
о 
(2-42) ee = t) 2 e E 
Let A - I 
2 
ала: 
DET 
r S - 
(2-45) Y= C z(Art) e 
where: 
il 
| 
E 
and Craig shows, that: 
| 2 
A 
(2-4) سف ے ن‎ 
eB 
e TA 


pem@ec iG 125 the purpose of this thesis to utilize 
Pearson's frequency curves in fitting the field data under 
consideration, rather than to develop the mathematical 
theory upon which the curves are based, which is a complete 
thesis in itself, only limited discussion of the development 
of the curves has been presented. For a more complete 
coverage of the development of these curves the following 
sources may be consulted. The foregoing discussion and 
derivations represent a synopsis of information given in 


the following references: 


a. Annals of Mathematical Statistics, Volume VII, 19356, "A 
New Exposition and Chart for the Pearson System of 
Frequency Curves" by Cecil C. Craig. 


DN ccuenewbuPyes end Corprelat2zon, by W. Er Eildervon, C. 
Cambridge, 1958, 
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c. ‘Karl Pearson's System of Generalized Frequency Curves’, 
by Arnold M. Wedel, Thesis, University of Kansas, 1948. 
d. Handbook of Mathematical Statistics, "Frequency Curves", 


by H. C. Carver, op. 92-119, 1924. 


The entire family of curves is shown in Table III. The 
main types are shown followed in order by the transitional 
types associated with a particular main type. This study 
has been limited to the use of the three main types, the 
transitional Type III, and the Normal Curve. These five 


curves adequately describe the field data used for this study. 


(o, 8) Chart 


In the course of the preceeding discussion a set of 
conditions for the various types of functions has been 
established in terms of o and 6, parameters which — be 
readily calculated. The numerical values of these two 
parameters determine the Pearson curve appropriate to a 
particular distribution. The conditions for each type or 
curve are summarized in Table Iii. 

An (9°, 6) chart which gives visual presentation of 
these conditions and an automatic means for type identifica- 
SN TS relatively asy to construct. 

In addition to the lines, 6 = -1, 6 =-1/2, 6 =0, 6 


= 2/5, and = о, the Chart contains only the curve: 
2 


(2-15) Rs = 15 (6+2) 





THE CURVES, 8 
e 
peu] EQUATION CONDITIONS - l| REMARKS 


“фо | Limited range, skew, 
a <6, § 4-2 usually bell-snaped 
2 





I = u * |but may be U-ahaped 
| Ye C (keri) (a,-t) — (2+252 ¥ "баз дент шүк жой; 
| Е 
= 0 metrical,usually bell- 
EI | °, i yet shaped but U-shaped 
уље (ty? ¡ -1<9<0,9F 7"  |when € 1,8, Range 
E | is (==, з) e 
Unlimited range in 
lone direction, skew; 
III bell-shaped, but may 
be Jeshaped. Range is 
( =Á , «OO je 
КИ Jesnapeds Range рст 
VIII an. infinite ordinate at 
TE -or~ = " rn, to finite ordinate 
pis ^^i) 2425)o4" -V(2$)(255) fat TI 
As Ё Оо Range from t = ra to 
IX -2ж -1 <8<6 infinito ordinate at 
Ут € (^ ,-t] А “ny (+25) ¢ és. ¡Pue J-shapede 
Az FO је 
X Кз, M Jeshaped with range 
У = A," > 4/ (е1, co ) 9 
Jeshaped with range 
ta, 
оз o Unlimited range 
IV. -m te 7 
уг с [онај та“ КЕ FIR ы y ы. e skew, bell=-shared. 


з #0 Unlimited range in 
V ај „ze ó >o one direction, bell- 
¢ (7 T 2:2: 95 (5==) _|shaped, Бапве 18(-Г,с2) 
‚nlimiteä range, sym- 


* | í * = 

ош y= C(t Zu) E T metrical, bell-shaped. 

Nore КЭ: Unlimited range, sym- 

mal Сет ey = Š = o metrical,bell-shaped. 
to Unlimíted range in one 

VI s 2 direction, skew, bell- 

Ye Ce" (e um | É >o shaped but may be J= 
š o, Z > 4 $ / $ z shared, 


As to J-shnaped with finite 
ХІ — {з er Ó<S e vc ordinate E $ =г:, 
pis с5а [2 32). +, Z X t$ 2 «; |Капсе 13 (ra, nr 


TABLE III Summary of Pearson System of Frequency Curves 
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on which the points corresponding to the Type V function lie, 
and the curve: 


(2-46) (2 + 35) = 4(1 + 28)? (20 4 6) 


On which the points corresponding to the distribution 


ENunetrons of Туре VII, IX, X, and XI are found. €F 


Construction of Ca 6) Chart 


Point a =0, 6 =О satisfies the conditions for the 
normal curve and is спе starting point LOr Constructing The 
graph. The lines 6 = -1, 6 =- 1/2, 6 =0, 6 = 2/5, and c 
= О, are easily constructed. 


For the equation: 


(2-15) a, E 


For the equation: 


(2-46) (2436) а, y ao 





Note when 6 = - 2/3, the expression (2+36) = O, therefore 


the line for this equation approaches - 2/3 asymtopically. 
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Figure 2 presents visually the family of Pearson's 
curves in terms of 2, and 6. It is a simple matter to 
determine the type of curve that fits any set of data 
merely by determining the & and ó values for the data and 
then entering the chart with these values. 

Tne subscript, B, on the chart refers to bell-shaped 
curves, the subscript, J, refers to J-shaped curves, and 
the subscript, U, to U-shaped curves. 

The points for ô < -1, correspond to no frequency 
functions, they fall in the "Impossible Area." 

ке ol GCesienaved as Neverotypic chose members cl s 
system for which the eighth movement failed to exist. (In 
such a case the standard deviation of the fourth moment in 


samples would be infinite.) This area is established by 


DE x. 
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СНАРТЕВ III 
APPLICATION OF PEARSON CURVES 


РОСС оп 

The brief theoretical discussion of the Pearson 
frequency curves having been presented it is now useful to 
examine how these curves may be applied to the field data 
under consideration. 

From Table I there are eighteen sets of data available 
штит (оте fields. The Pearson system of frequency 
sI eS will be applied to find a good theoretical fit for 
each given observed distribution. 

Шшас 15 спе value of frequency curves? ~Aymormal curve 
fitted to a given set of data is to determine whether or 
not ое часа are normally distributed, ІТ the distribucion 
normal then use may be made of the extensive body of 
sampling theory applicable to normal populations.  Compara- 
tively little is known concerning sampling from non-normally 
distributed populations. 

When data are distinguished as non-normal, there may 
be further advantages in fitting a non-normal curve to the 
data. Such a curve may serve to smooth the histogram and may 
thus permit a more accurate determination of the relative 
Поа yC eS of che population from which the sample was 
caken. The identification of the distribution of the given 
data with a particular frequency curve may also serve to 
(ен Е бей (rom other data, the distribution OF ииси 


is identified with a different frequency curve. These are 





> 
the two principle reasons for fitting non-normal frequency 
EuPves. Also, if a particular type of non-normal distri- 
buvion oecurs with sufficient frequency, this fact will 
Serve as an incentive for the creation of the appropriate 


sampling theory. 


Fitting the Data 


wine arichmetical labor involved in fitting @ sev of 
observed data to a frequency curve is lengthy and tedious. 
This may be a prime reason why there has not been greater 
utilization of Theoretical curves in practice to obtain a 
description of the distribution of a given set of data. 
With the advent of hish speed digital computers this 
objection to the heavy arithmetical work involved is lessened. 

To illustrate the procedure involved in fitting a set 
of data to a Pearson frequency curve, the step by step 
calculations for one set of data, that of permeability for 
Field 3, is shown. The calculations for fitting the remain- 
ing sets of data were performed on a digital computer. 
The computer programs, written in Fortran language, with 
accompanying Flow charts are shown in Appendix B. For each 
data set, the histogram and plotted curve follows the 
numerical calculations for that set. The normal curve has 
also been fitted to each set of data to illustrate the 
comparison with the non-normal curve that describes the 
data. In many cases the goodness of fit test for the normal 


eurve indicates that the normal curve gives a poor fit for 
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tne data, which is as one would expect for some of the more 


markedly skewed distributions. In a few cases though, the 


normal curve does give a satisfactory fit to the data, 


even though the normal fit is not as good as the frequency 


curve describing the data. 


Procedure 


The procedure used in fitting the various types of 


curves is: 


n. 


“м 


Ол 


95 
Ls 


Arrange the data in an array. Tabulate the data using 
convenient class intervals. 

Calculate the first four moments about a convenient 
vertical. 

Transfer the moments to the centroid vertical or 
vertical through the mean. 


1 to the moments if there 


Apply Sheppard's corrections 
is high contact at both ends of the curve. 

Calculate E фу, and o. 

Ic rese mean X and the mode U. 

Determine by use of the C о) саг лш но 
curve, Should be used. 

Calculate the constants for the equation ог the curve. 
Calculate the theoretical frequencies at the mid-point 
of each class interval. 

Plot the histogram. 


Plot the theoretical curve constructing the mid-ordinates 


at the middle of each class interval. 
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l2. Calculate the area graduation for each interval. Test 
ror goodness of fit of theory to observation. 
15. Fit the normal curve to the data for comparison and 


test for goodness of fit of the normal curve. 


sample Calculations | 

The permeability data for Field 5 consists of 195 
observations obtained from seven core analyses from seven 
wells. These data were divided into 14 class intervals of 
10 millidarcies each. Table IV shows a convenient method 
of calculating the first four moments for the core data. 
Following Table IV are the calculations needed to fit the 
frequency curve to the data. 

The fitting of a set of data to a Pearson curve is 


th 


based on the method of moments where the r moment around 


anar bic rary origin 15:2 


pom) utc Г х’ ав(х) 


—со 


The zero-th moment about the origin, Ms» always exists and 


is equal to one. 


The pen moment around the mean is: 


(e) Ln = INE. uj!) ar(x) 


009 


If the moments around an arbitrary origin are known, 


the following formula is used to find the moments about the 


2 


mean: 
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Б -2) Еи ЕЕЕ | 


"us: 


2 
Ha = Ho’ Mi’ 


= қ > 
Hz = Mz | لار‎ H. ' + Cp ' 


г | д 
Hy = Ly | = Hun Bs ' + бш 14 B. ' = SH! 


се ве ol and calculation and simplicity ornnotrarion; 
the scale for the independent variable of the frequency 
Bsopibpubtjon, (i.e6., porosity, permeability, or saturation) 
was transformed to a notation wherein the interval contain- 
ing the arbitrarily chosen mid-point is numbered zero and 
intervals on either side are numbered serially. The 
negative values were assigned to that side of the distribu- 
tion wnich contained the mode. This convention is in keeping 
with that adopted for the sign of "a" in the derivation of 
the Pearson system of frequency curves. An additional 
transformation of the independent variable is made to correspond 
to standard statistical notation. This transformation is 
the introduction of the standard unit t. Thus when the 
histogram intervals are numbered serially, 


P 
2 L 








En 
where x, takes on the values 1, 2, 3, 4, etc. at interval 
mid-points. 

The equation relating t and the original physical 


Everly, d, is: 
(3-5) а= $ о да+а 


where q is the mean value of q and Aq is the interval of q 
used to construct the histogram on the q Scale. 

ine following calculations starting with Table IY are 
for the permeability observations for Field 3. For each 
of the other data sets only the resulting constants and 
equations are shown followed by the graph of the frequency 
eve and normal curve. 

Taking moments about the arbitrary mid-ordinate 45.0 
in terms of the transformed variable x with its correspond- 


ing mid-point, x = О, the following moments are calculated: 


А EM И 
(5-6) из! = ши“ = 195 = -. 46152846 
(3-7) | EX 21260 6 1615585 
27 a LOS 
(3-8) 227 960 


м5! = = = == = 131. 2000 





DO 
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{86 ‘ве 096 - 091 = 06- - GGT - E MA 
199 ‘9 1999 627 Ж 755 Fo 6 I GET Oy T=-TO*O€T 
960 *+ 0604  <Tú ZIG +9 "9 9 9 І GZT O I To O í 
ток с ТОК СС HE бї Өк m D t атт Os (LOMO 
обет обет 912 оте 9¢ 9€ 9 9 [i | GOT Sao oc 
OGZ ‘T сад 092 чет OG Go OT G 2 аб | OOL= TO 06 
3403 оба ата 119 Зет e 7 1 8 48 06-10 08 
168 t9 169 6c 66 6 сс C ET. . GJ, Doe 
ода OT бет 9 9 H po а 9T G9 TOUS 
ог fi O2 t ог її O2 t O2 GG 09 -I10'0G 
О 0 0 0 O O О О пе SH OG -TO'O 
zç Í oC - T- гс il Zg- T- Zg Gc Оң -І0704 
күб от слан OS “аст H 00- Zz- КС Ge 06 -то'ог 
CHEZ 18 €9l- L2- тоё 6 Lg- g- бе GT 02 -TO*OT 
048 “E оба 096- 9- Ote ao mR GT G ОТ Б 
хх A J | A Jx „х ух x E (О тл; jo STADE 


Kouonboua, эптел-рти Азтттаеэшлэа 








© АТЯТя--АТТИТаАхнична JO NOLLOSINLSIG WHL ЧО SXNHHON NOI LSUTA HHL HO NOLLVINITVO 
AIL WISV | 





39 


And, for moments about the mean, 
Ha = 0 
lo =p! ш! = 6. 1615585 - (-. 16155856) 


Ha = 6. 24185208 


1 I t j > 
Hz = Н» = ÍH! Ho w e» 
из = 1.9220769 - 3(-. 16153816) (6. 461538) + 2(-. 46153846)? 


2 4 
My Ay! - 4! ы. ou! Ho! = Ju! 


|| 


Му 148. 41116 


And the standard deviation for grouped data, 


(3-9) о =Уй> =/5.1651875 = 2.1997055 


Hence Sheppard's corrections are: 


(5-10) X, = Hu - == - 6,2485208 - .085555 
My OESTE 
(5-11) X, = B, = 13.673191 
y NE 
El Ху = № - 3 Mo * 250 


X, = 148.41116 - =(6. 2485208) + .029167 


Ap = 145,31607 


us 








10 





Ку 14551607 145.31607 И 
фу = x2 TS = 3.8231470 
5 2 = - 5 а, - 6 
o ЕД + 5 


5 . 2(3.8231476) -3(.79781261) - 6 
(3.8231476) + 3 


ö = -.10950116 
DLS 2, - 45(6+2) 


= „79781251 -4(-.10950116 ) (1. 8904988 ) 
= 1.6258600 


Then entering the (о, 6) chart with: 


a= = .79781261 and 6 = -.10950116 


р 


it is observed that the frequency distribution of the 
permeability for Field 5 may be fitted by а Pearson Type 
lg curve. 

The equation of the curve is: 


m 


m 
1 (2-6) 2 


à C(t-r, ) 


where: 
N Г (ту + м. + 2) 
С = 
И 
= r) T(m +m, ) P(m,+1) 


where N is the sum of the frequency for all class intervals 


and: 


Е + /D 


-— Aas 
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p. 27299320367 4 1.625350 
s 2(-.10950116) 


vr], = -1. 7437833 


r... 
еее 
2 E 


г. = 9. 9007883 


s. s аш 


JD 
NE 900950 DIS . 89320367 
ту = - So) - 179976919) - 1 
ту = 1.4356194 


| а 
m, = EEG + 2 -1 


m, = 12.829027 


And: 
C = 195 Tr(1.14556194) + 12.829027 + 2 


(9.9007883 % 1.7157655):2:29%955 11) 264646 г15.829027 


195 T16. 2046464 


az 15.29 
(11.64415516)^-* 113 7 56161. 29027 
алтоо ри 
The Gamma function, defined aes 
(5-12) r(x) = f a 60-1) ас KG 
О 


was calculated in two ways, depending upon the value of x. 





12 
For values of x less than 9.0, the Gamma function was obtained 


from the recursion formula: 
(5-15) ГО np o) 


which is useful for expressing T(x) as a function of some 
value of T(a) where (a) is a value between 1. and 2. The 
value of T(a) was. determined by a table look-up and 
interpolation using tabulated Gamma ee 

For values of x above 9.0, P(x) Was computed from the 


T 


Soreling formula’ for log r(x). 


log T(x) = -.15429448(x) + (x -.5) 10510 X 
+ 108,0 (1 + $22] + .39908993 
леса оп ог tne curve is: 


-12) l. 4356194 12029 


Ш- (2.211844Х10 (t+1. 7437533) ПОГОН) 


The range of the curve can be computed from the above 
Ес Lion. 

The range of the. curve is defined as the upper and 
lower values of the independent variable between which all 
Positive frequencies exist. These theoretical limits are 
determined by seeking the values of t for which Y = 0. Ву 
inspection of the above equation, Y is zero when t is -1. 74576335 
and +9.90070. These values of t correspond to permeabilities 
of -2.65 md and 288.35 respectively as determined. from the 
equation. 


q = q + t сада 





25 
where t is one of the limiting values given above, o the 
us rd deviation in units of x, А q is the elass interval 


165 of q, and q 


= 


is the mean value of q. Because negative 
values. of permeability are not physically meaningful, the 
limits can be looked upon as zero and 208.35 millidarcys. 
Ine existence or the negative lower limit should present 
no more interpretational difficulties than those presented 
by the known limits of any normal frequency curve which are 
by definition, -o and +, 
For calculation of graduation (mid-ordinates), the 


rollowing arrangement is convenient: 
(1) (2) 


Distance from 
Origin 






Mid-poin 





2,6153846 . 984173219 |2. 7284955 |. 15592 18. 91605951 


(т) | (8) (9) | (10) i) (12) 


Area ог 
Тог (6) m,x(5) mox (T) (3)+(9)+LogC=LogY, Y. |interval 


ы. 


‚ 95017 1.624307 | 12.18975 | 1.16038 | 11.167] 14. 512 


The area of the interval is found by means of Simpson's 
во 


ү: 
(5-15) f £(x)dx = 1/6(Y, + ^Y. yo ae 
O | 


By use of the above calculating techniques, the values 
of Y,. (Graduation-mid-ordinates) and Graduation (Areas) was 


found for each interval. Additionally the values of i OI 





11 
(Graduation-mid-ordinates) was found for each mid-point by 


substituting the appropriate value of t for each mid-point 





Шо the equation: 5 52 
» P 
SM ON e 
where: 
on ЕД N 
мот G 
Thus: E | 9847 2 
О 
Yan = 2d e 2 Сл 
> 
гт О 


The following summary for Field 3 permeability data is 


presented: Type 18 Curve 


ү = (2.21184 х 10712) ($ + 1.7457655)1- 256194 (9.90079-5)12.829 
Permeabilit I Graduation Graduation Normal Curve 
fees int (may ET (mid-ordinates) (Areas) (mid-ordinates) 

5 15 14.73 14.48 11.43 

115 29 29.14 23.72 gu 

25 34 33.85 33.57 25.75 

p EA Due Sloe 30.41 

1.5 24 26.47 26.44 30.60 

55 20 20.23 20.25 20.28 

65 18 14.47 14.51 19.16 

75 11 9.77 9.82 c 

85 3 сас G ОЕ 

95 2 3.82 3.86 2а 5 

105 1 2, 22 2.24 12:10 

Du 1 1.22 1.24 NO 

125 1 . 61 ‚65 TIO 

135 db . 2 02 

195 194.95 194.09 184.36 
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Goodness of Fit 


The Chi Square (72) goodness of fit test for the above 


curve may be calculated by the following method: 


Chi-Square "Goodness of Fit" 





P 1 f- F (f - ғ)? Im 5 D 
15 14.5 H5 8-5 2015 
29 2957 22 .09 .005 
p 5580 E . 16 . 005 
в = p 53 . 09 . 003 
24 26.1 -2,4 5.76 . 225 
20 ВОДЕ - .5 199 . 005 
16 UNE 1.5 2.25 .152 
ial 9.8 1862 1.44 ‚142 
8 51% iem 2.89 ¿444 
2 Бо -1.9 3,61 . 895 
1 2.2 -1.2 1.44 .689 
d D - .2 . 04 . 046 
al ҰЛ p . 09 105 
Jj iE E Sot . 700 
194.1 2, ӘРЕ 


The number of degrees of freedom for this curve isn - 
1 = 15, where n is the number of class intervals. Entering 
a X^ table? with 13 degrees of freedom and X? - 3.522 an 


approximate value of p = .99 13 found. This indicates 

that there are 99 chances out of 100 that differences as 
large as those found could have arisen due to chance or 
Sampling variation. Therefore, we can conclude that the 
Type La gives an exceptionally good fit to the permeability 
data. For testing the goodness of fit the significance 


level of p is usually taken as either .05 or „01.29 
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Numerical Results 


Permeability (K)--Field 1 


Data Range (.2, 439 md.) Mean = 63.26 md. 
Class interval = 40 md. с = 60.8 па. 
EC = 3,0327404 5 = -. 2031506 Type I, curve 
г) = -. 95047725 г. = 9.1901735 
m) = -. 3008110 m, = 6.1312036 
= 2.8551007 x Io 
Y = (2.8551007x10”*") (t+. 93047735) 730081100 ¿(ç 4901735-t)9- 1512056 
Curve Range (6.7, 490 md.) 
15 Curve Normal Curve 
ат X^ - 968.2 
ENS good at .O5 level Pig 16 ПОШ 5006 
Permeability Graduation Graduation Normal Curve 
E Frequency (mid-ordinate) (Areas) (mid-ordinate) 
20 ЩЙ 610 786.1 600.8 253.3 
60 281 Ў 520.9 329.7 - 558.7 
100 189 174.1 O 291.6 
140 101 97.7 99.1 106 Y 
180 a 54,4 55.1 SUS 
220 25 29.4 29,9 14,7 
260 ШЕ 15:2 m 22.0 
200 9 T5 do . 2 
340 1 3.4 3.5 . 02 
380 E 1.4 1.4 0007 
120 1 „5 „5 00002 
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Porosity (¢)--Field 1 


Data Range (5.4, 25.5%) Mean = 19.45% 
Class interval = 2.5% б = 3.65 
a, = „72488482 ö = -, 29374565 Type In curve 
r, = -1.3630479 Го = 4.2614796 
my = .16551900 m, - 2.6432920 
С = 5.7594065 
У = (5. 7594063) (&+1.5650179)1.6221900 (1, 0614796=t) 29452928 
Curve Range (3.83, 24.0%) 
18 Curve Normal Curve 
X* - 19.841246 Xe ede 
Mil is good at .O2 level [spen M PESE ODE 
озу Graduation Graduation Normal Curve 
a оо (mid-ordinate) (Areas) (mid-ordinate) 
6.25 + . 220 2T . Д 
8.75 5 10.9 11.4 DRI 
11# 25 32 DIS DS HO © 
15. 76 34 64.0 64.5 67.4 
16.25 90 EON TO G 15140 
y 159 Loom 168.2 213.6 
21.25 278 228.8 Bo See. 
Ze 140 246.9 206.6 104.6 
26.25 5 on ој 36.4 


794 
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Oil Saturation Е 1 


Data Range (7.0, 73.0%) Mean = 38.89% 


Class interval = 5% б = 11.05 
E = .020572560 5 = .027215451 Normal Type curve 
С = 142.75 
ce 
ШЕ 1125 75 - 2 


Curve Range (-e, e) 
Normal 
Ye = 16.568630 
lues rood at .25 level 


Кыс = 
О 


Normal 


Midpoint (%) — NOY Mid-ordinates 
8.5 6 UNA 
15.5 20 11.5 
18.5 29 28.3 
25.5 48 57.7 
28.5 ЗД 95.8 
22005 116 129.7 
38.5 162 143.2 
13,5 141 TOS 
43.5 98 94,5 
55.5 55 56.5 
58.5 20 DO 
65.5 12 11.0 
68.5 1 Ет 
ов Et . 2 
794 
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Water Saturation (S,)--Field 1 22 


Data Range (15, 89%) Mean = 40.96% 

lass interval = 5% о = 12.15 

а, 2725375019 6 = -.10556269 Туре I, curve 
ry = -1.8241159 E 10.038806 

m, = 1.6619956 m, = 13.649978 


С = 9.5915501 x 10712 
y = (9.5915501x10"5) (t+1.8241159)1:661995 (10, 938806-4)13- 249978 


Curve Range (1850708 


In Curve Normal Curve 
y? = 10.875972 X? = 95.56815 
EC 15 good at .80 level HI 1S no wege 


W 22 Graduation Graduation Normal Curve 
Midpoint (%) ~ У wid-ordinates (Areas)  Mid-ordinates 


i 8 О б 1246 
21.5 27 35.4 36.3 38.6 
26.5 39 = MS ES 
Е Ш.Б 150 136.5 d il 99.2 
20.5 160 136.8 156.0 123.4 
41.5 136 Du s 117.5 jc 
16,5 78 91.6 ee 115:2 
51.5 51 65.9 DO 86.3 
56.5 ДО 44.3 44.6 54.6 
ine 52 20.1 po 29.2 
66.5 15 15% IO 2 
Dres 10 9.6 ү 5.0 
Hons 6 s 5.2 MS 
81.5 5 26 В.Т .Д 
96.5 2 12 1605, zu 
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Perneability (K)--Field 2 


Data Range (.16, 381.0 md.) Mean 


И. 5509951 x tone 
Y = (1.8509951x10”*) (t+. 71192903) ” 
Curve Range er 
lj Curve 
Ж = 21.96777 
Fit is good at .02 level 


= 49,5 md. 
Class interval = 30 md. y ПО. 
| I = 5.4370534 ö = -. 22554678 Type I, curve 
гу = -. 71192903 Po = 11.050699 
ту = -. 58435670 m, = 5. 4516905 


. 5343 5670 (11.050699-6)2° 516955 
636.8 md.) 

Normal Curve 

у“ = 486.197 


Ric 182100.2063 


[| ү; Freqieney md ordinates (Areas) (та NE 
15 295 1080.9 Э AS 
1.5 164 SIO 159.1 О 
75 ҮТ 66.8 O 1 1284: 
105 52 39.4 39.8 ОБЕ 
155 28 ale 24.6 41.8 
165 12 HM 15.5 1 
195 iy Sic di 9.8 Dac 
225 5 en E E 
255 2 D e : 
285 n 2,2 2.5 . 01 
ENS О D > .0007 
DE 2 T Ti . 0000! 
375 1 4 4 .000001 
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Porosity (¢)--Field 2 


Maca Range (8.1, 24.3) Mean = 19.72% 
Class interval = 2% о = 2.34% 
а. = 1.3035387 6 = „14055615 Type Vi, curve 
шу = 4,5488532 m, = -22.778043 
E- -1. 9595330 ro = -7. 7511888 
С = 20.968745 
Y - (20.968755) (&&-T. 7211868) 29: T1605. (419698380) 979522 
Curve Range (-o, 24,33) 
Via Curve Normal Curve 
X^ 2 15.456512 X* = 71.647549 
Fit is good at .05 level Fit is not good 


Graduation 


EOrOSLCY 5 
Y (mid-ordinates) 


Midpoint (%) *Yequence 


9 2 vi 
11 2 1.9 
13 6 5.4 
15 9 15.0 
"a 47 HO. 4 
19 106 9h. 4 
2 159 ШІН 752 
25 38 DR 
25 E 0 


N 
I 
JI 


Graduation Normal Curve 


Areas (mid-ordinates) 
"M. 20057 

270 sd 
5.6 1.9 
и 15.4 
11.7 59.6 
94.8 ae 
AS TOL 
17,0 2.1.5 
2 9.5 
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О Бане о оп о 2 


Data Range (9, 79%) Mean = 49.38% 

Mass interval = 7% | o = 11.7 

а = .6911813 6 = -,30505016 Type I, curve 
г] = -1.2600192 Гр = 4.0852955 

m, = .15792540 m, = 2.1181582 


С = 1.749383 
у = (1. 7198815 ) (5+1.5600195)1-27925* (1.0852955-5)2- 1191592 


Curve Range (-10.7, 69.4) 
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Із Curve Normal Curve 
x = 21.757241 m = 90.32 
Fit is good at .Ol level HIC iS OL оо 
S ` Е | 
Graduation Graduation Normal 
Midpoint (%) Frequency wa opdinate Areas (mid-ordinates) 
йү 15 ои! 6.1 2.6 
19.5 T2 11.4 qu O 
25.5 5 119156 ПО | йо с 
ре. 5 21 205 IR 55 
99.5 pal por DO Doe 5 
15.5 Sal 50.5 50.6 64.7 
53.5 71. 635.1 O 2255 
60.5 76 74.7 74.5 17.9 
[Б 29 y occ 61.4 29.4 
74.5 2 O. e 14.5 
Q1.5 l 0. О, 5.6 
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Water Saturation Пе ено 2 


Data Range (11, 85%) Mean = 38.65% 

Class interval = 6% . c 214.82, 

E = 1.0579891 6 = -.32958992 Type I- curve 

ШІ - -1.178(589 А2955 Би 

my = -. 12466560 m, - 2.192812! 

C = 1.8860123 

Y = (1.8860125) (6+1.1787589)7: 12:00500 (1 29956344) 2-192812% 


Curve Range (21.2, 102.3) 





15 Curve | Normal Curve 
X^ - T.3586261 X* - 98.996255 
о је good at .75 level FILES HOL 2O00 
о UPS Graduation Graduation Normal 
Mid-point (4) амелеу wid-ordinates Areas (mid-ordinates) 
2, 1 Qu O. J z 
19 18 E D 25. 2 
25 75 SES 82.9 D D 
ei 12 О oles 19.5 
E Sl 5 47.6 ا‎ 
34 36.9 p 52.8 
19 25 207 28.4 12,6 
ES 18 21.2 21.9 29.2 
61 16 15.4 1545 17.09 
67 9 1047 LO 7 8.4 
12 2 7.9 7.0 3.5 
19 6 4.1 4.2 L25 
85 E 221 a I4 
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Porosity (%)--Field 3 


Data Range (10.2, 23.1%) Mean = 18.06% 

Class interval = 1.5% о = 2,46 

a = 581951115 ö = -. 25091835 Type I, curve 
p. - -1.526424 р. = 4.5666935 

m, = . 4957617 m, = 3.4749591 

Ви . 11555541 

Y = (.11555541) (t+1. 526424) "+ 957617 (1 5666935-t)2: "79591 


Curve Ran bv ЗА и 


lp Curve Normal Curve 

e 752871 X^ - 36.031158 

в 15 good av .O5 level Fico БКО ОО 
woROosi uy IO ERO Graduation Graduation Normal 

Midpoint (%) equentY mid-ordinates Areas (mid-ordinates) 

O 75 2 ШБ Ш . Д 
12.25 2 До Sat t 
d (^5 ae ТО ТЖЕ 1 
15.25 l2 1666 1522 1720 
IT 20 UNS Вахо 28.0 
2.25 | 25 29.0 28.9 31:9 
LSO 48 о ES au! 
ale 5 9 25.5 22.5 19.0 
22.75 3 O O. E 


192 
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Oil ласты ко Ж Есе Б 


Data Range (13.1, 48.7%) Mean = 33.26% 


Class interval = 3% © = 5.7 

ЖЫ - „79434377 6 = .O72484741 Type VI, curve 

ру = -1, 2828677 рр = -6. 6759047 

m = 51. 961440 m, = -83 . 553454 

С = 37. 514620 

Y = (37. 514620) (t+6.6759047)783-5555% (64 0828677) 51: 0 


Curve Range (-», 57.7) 
Via Curve Normal Curve 
Ж - 10.861505 Х2 = 24.761836 


ІІ 15 соса at .50 level Fit is good at .Ol1 level 


5 


Midpoint (4) #requency TT e a 
14.5 Í: ES 7 al 
Jm S Y 5 1.4 26 
20.5 2 2.9 2.9 2.4 
22.5 2 5.8 DS One 
26.5 T: D a 1120 
29.5 18 19-7 те 22.6 
32.5 210) 26 15 28.6 27.8 
pou 54 209 29.7 D DEEP 
EC 19 2211 DOM erde 
41.5 9 Je 9.9 9.4 
14.5 | 1.5 1.8 n 
17.5 1 04 08 ШІ 
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6T 
Water Saturation пе © 


Data Range (30.4, 76.1) Mean = 45.35% 

Class interval = 5% o = 7.85 

a, ^ = 1 46519% ô = -.039926Ł2 . Type III, curve 
A = 1.6534022 Е. 

Е СОА. QA 71 Si 

Y = (13.616174) (1.6554022)7-38950 «71.0534022 t 


Curve Range (32.4, c) 


Ill Curve Normal Curve 
X° = 4.7735815 X^ - 16.889006 
Bit is good at .85 level Fit is good at .05 level 
W Graduation Graduation Normal 
Midpoint(Z) Fremencey wig-ordinates Areas Mid-ordinates 
Do. S 6 qm 4.1 9.3 
p. 20 E BBS ey 
42.5 13 ү 2 Тб EIS 
17.5 29 205 26.3 320 
E25 He 14.9 I 21.4 
E Ө p T 9.6 
ОО», Y 3.6 sn 2 
oS 2 06 nd 15 
2.5 О 27 M ET 
55 1 5 Bo 
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69 
Porosity (@)--Field 4 


аса Range (9.5, 22.4%) Mean = 17.8% 

Class interval = 1.5% С = 2.0 

a, = 1, 1315972 ô = ,0047189 Type III, curve 

A = 1.8801133 С = 18.084424 

2 4 
о) V 
Y = (18.684424) (1.88011334t)2 25440 ¿-1.600113 t 
Curve Range (-», 21.57) 

tit. Curve Normal Curve 

X? = 6.303058 IS 

Fit is good at .60 а 

Boros. UY | Graduation Graduation Normal 

Midvoint (%) Hrequency wig-ordinates Areas Mid-ordinates 
9.75 2 2% . 9 05 
MaS 2 2.6 2.7 ED 
вт 8 атар Le 4.0 
14.25 16 os © Im 20 
15.75 36 41.3 11.9 57.0 
Е 25 95 79.0 Tone 91.7 
19.75 105 108.2 104,9 03.9 
20.25 53 ОЕ 59.2 ЗЕ 
21.75 5 О as ША 5 
2050265 О О OR 232 
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SETI о П 


Maca Range (15.5, 68,2%) Mean = 59.96% 

Class interval = 3.5% o = 10.1 

2% = ,01836461 ô = -,13178302 Type L, curve 
ру = -3. 2859150 a 

пу = 4. 6958277 m, - 6.796547 


1.1959525 х 1072 


|| 


Y = (1.1959325x10">) (+3. 265945) *+ 298277 (1 3110704-6)0- 4796547 
Curve Range (6.9, 83.3) | 
ја Curve Normal Curve 
Y^ - 13.992916 X* = 15, 734882 
Fit is good at .60 level Fit is od acer EN = 


Midpoint %) ... Б С и ail ОН 
14.75 n . 99 N IC 
15.25 6 Ds 4.0 1.1 
21.75 2 9.6 9.7 ЭЎ 
25, 25 au DS Г ШЕТ 
28.75 21, 26.6 26.6 25.9 
Be. 25 ES о Pi 34.5 pone 
55.75 15 10.0 59.9 10.5 
20.25 39 41,6 11,5 15,8 
12.75 32 DOS ES 42,1 
16.25 29 БЕН 2259 Bion 
19.75 25 26.5 26.5 PD 
55.25 26 о 1007 1340 
56.75 5 qr DNO 10.6 
60.25 T 6.4 6.5 5.6 
03.75 3 510 eo ^ Ne 
87.25 2 т ze па: 
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Water Saturation (S )--Fleld 4 


Data Range (18.6, 76.84) Mean = 59.81% 
Class interval = 6% б = 15.5 
a = ,44508734 | 6 = -.41374306 Type I, curve 
ry = -1.3112949 р. = 2.9237671 
EN -. 1202559 m. = . 956141568 
= 10.194623 
Y = (10.194628) (t+1.3112949)7: 12259 (2.9257671-5): 9501568 


Curve Range (22.1, 79.2) 


ig Curve Normal Curve 

У“ = 15.106190 YE = 28.03312 
ШО iS соса at .05 level ULIS E CU DC 

W ee, Graduation Graduation Normal 
Midpoint (%) Brequencey wig-ordinates Areas Mid-ordinates 

21 19 O 220 14.2 

2'{ ДО 41.3 41.7 25.4 

22 45 2222 2988 puse 

1.5 18 250 0E 320 

51 20 Ше 5 HM 21,6 

E 20 14,4 14.4 2532 

65 T 10. 4 10.4 К 

69 б 6.5 a e 

75 3 2.6 236 jet 

201 
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Porosity (@)--Field 5 


Data Range (3.4, 34.3%) Mean = 20.66% 

Class interval = 3% б = 41,32 

Же = „33849171 ö = . 25950502 Type IV curve 
р - 1.1205492 m = 5.8520056 


У = -5.9842952 2. (291120 


|| 


|| 


33495724 х 10! 


У = (5.510572х107) [ ($41.1205892)2 + (2.7291720)2]-2:8520055 x 


n (StL. 1205792) 


го веат "n 
3.9042952 tan 2. 7291720 


° 


Curve Range (-o, о) 


IV Curve Normal Curve 
о > u 
x = 83.10963 Ne = 109. 93954 
ШЕ 15 пос good FiG 15 nov ова 
[IODOSlcy ara: Graduation Graduation Normal 
Midpoint (%) Guen’Y wid-ordinates Areas Mid-ordinates 
4,5 9 1.2 4.5 4.5 
E 155 os laka 1.7 
5 59 25.5 E 205 
125 96 99.5 10590 T2079 
шо TOn 245.4 2418.6 296.2 
1 472 150.6 143.5 249.0 
pon 592 — 170.3 120.3 
25.5 192 257.0 259.9 242.9 
20.5 LO о Г 36.8 
pics 29 Bor er pow 
34,5 1 n ee 2.0 
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011 Saturation | eld 5 


Data Range (5.7, 77.2%) 


Class interval = 8% 


E" = ‚00047316776 
E = -3 . 9580728 
m = 6. 5279175 


СИЕ 5711286 x 1079 


на 
|| 


Curve Range 


X^ = 16.2566686 


Fic is good at .05 level 


So 


(2.5711286х1079)(:-5.9500728)0- 5279175 (1 1474939-8) 


Graduation Graduation 


Mean = 32.77% 


o 15.04 

6 = -. 11483628 Type I, curve 
Гр = 4. 1471939 
m, = 6.8881809 


5.8881809 


(=21.9 90.1) 
Normal Curve 
X= = 17.879294 
Fit is good at .0% level 


Normal 


Midpoint (%) Brequency yig_ordinates Areas  Mid-ordinates 

д 65 TONO 50.9 16.1 
2 178 141.0 uo 128.5 
20 189 263.4 262.0 256.1 
23 343 555.6 549.9 366.0 
36 399 pons DE CO poo 
En 297 27586 272.1 2 
52 156 153.7 154.8 142.8 
60 48 58. L 60.6 55.4 
68 14 12.5 13.9 14,3 
75 5 . 9 D 2-й 








TO 


d 





тешлом — — — 
Т ЭЧАТ ——- 


NOLLVANLVS TIO 


% OS 


(7:06 “67182-) Омун 


Er TEAL OL AUOT 


9| 8 


001 


002 


00€ 


ООР 


Aduanbal4 





се балт OS Ес 5 


Mate Range (12.4, 92.35) Mean 52.17 
lass interval 8% o = 10.0% 

2 
a" = ¿000611895 ö = -, 02569028 Normal Type curve 
С = 244,82 

+= 

Ben} a 

Y = 244,82 е 


Curve Range (-, e) 
Normal 
Х2 = 14.359695 


Шис 15 £000 at .20 level 


о АЕ "m Normal 
Midpoint (%) pue ie Mid-ordinates 
12 155 S 
20 10 EOM 
20 56 72.7 
26 111 156.3 
Ts 257 198.8 
52 oo — 
60 189 | 199.2 
69 106 TOS 
75 76 5.2 
E д 30.4 
92 10 EI 


11529 





Water Saturation в ет 5 


Data Rance (12.4, 92,50) Mean 52.17 
Class interval 84 o = 16.0% 
a - :000611895 8 - -.02569028 Normal Type curve 
С = 241,82 
EU 
Eo) 2 
Y == оц Во e 


Curve Range (-e, о) 
Normal 
2 P 
X = 11.559695 


Fic is good at .20 level 


S 


wW - Norma. 
Midpoint (%) TT Mid-ordinates 
o 15 9.7 
20 10 SOR 
20 55 Саи 
36 ji 136.3 
11 257 198.8 
52 21:2 225.6 
50 189 | 199.2 
69 106 136.9 
76 76 oe 
gu 17 30,4 
92 LO 9.9 
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Summa 

A concise and easily visualized summary of the distri- 
putions of the data may be shown by means of the I 5) 
chart. Figure 2 is the complete (9°, ô) chart, бош бен, 
for purposes of showing the distributions of the properties 
under consideration, only the necessary portion of the chart 
has been prepared for each of the field properties. 

Pigure 21 shows the distributions for porosity for the 
five fields. There appears to be no pattern to the types 
Пс зрело јорз, іп that Fields 1 апа 5 are of the Type Ip: 


Field 4 is a Type III 


В? Field 2 is a Type Vi, and Field 5 


15 а Type IV. 


>< 
3 
0 | 2 
1; 
Е 
° 2-3 
— .2 
Те 
5 em ,] 
о д Ша 
2-4 
4 | Mg 
= .2-2 
.2 Na 
° 8-5 


3 


FIGURE 21 (a,*; ó) Chart of Porosity (Normal Curve is Pt. 0.0) 


e^ 





82 
Of the five fields, only Field 4 satisfied the goodness 
SI TIL test for the normal curve. This is somewhat in 


Hun 


contrast to the conclusion of Jan Law that: 


"With some exceptions permeability and porosity 
assemblages give respectively satisfactory logarithmic 
and arithmetic normal frequency distributions." 

For the present study it was the exception where 
porosity data gave a satisfactory arithmetic normal frequency 
distribution. . 

The (o, *, 8) chart for the S, data, Figure 22, indicates 
a somewhat different situation regarding normal distributions. ' 
four of che five fields, 1, 3, 4, 5, satisfied the goodness 


ЕП fit test for normal distributions. 


exc 4 
о г 
— а 
So-2 1, 
= ,3 е 
— .2 
50-4 La 
9 – „рео-5 
о š Ils 
So-| 30-3 
Ix 


FIGURE 22 ^, 5) Chart of S, 
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All five fields for S were fitted by a Type I curve 


2 


with varying values of o and 6. Field 5 satisfied the 


goodness of fit test for the normal curve. 





FIGURE 23 (с, 5) Chart of S, 


It is noted that both the So and Si distributions for 
Field 5 which has a considerably greater number of samples 
included, approximately. 1700, than any of the other fields 


satisfies the goodness of fit for the normal curve. This 


could indicate that if large enough samples are taken, the 


So and Sy distributions approach normality. For sampling 
considerations, the assumption that бо and S. populations 
are normally distributed would permit many established 
питон LO be applied го the analysis of these properties. 
A study was made of what effect the choice of interval 
size would have on the a and ô values. | For each set of 


data various interval widths were used, i.e., for porosity 
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Шрегсепс, 1.5 percent, 2 percent, апа 3 percent interval 
widths, to calculate the parameters а, and ó. For all 
fields with the exception of Field 5, no appreciable dif- 
ference was noted in the values of i and ó for the various 
interval sizes tried. There were only 132 observations for 
Field 3, and the calculated values of a, and ó were more 
sensitive to the size of interval used, with а, ranging 
from .641 to .279 and 6 ranging from -.249 to -.099 as the 
interval size was varied from 1 to 2.5 percent in steps of 

. 5 percent. 

lt was concluded from this study that with a sufficiently 
large number of observations (> 200), a width chosen to give 
between 7 to 15 intervals will give consistent results in the 
calculation of Ong and ó. 

Although it was not possible to study enough fields to 
determine the relation between reservoir type and the 
parameters а, and ó herein | it would be reasonable 
to assume that such a situation exists. Specifically the 
parameters а, and 6 would be expected to reflect the 
variation of individual properties throughout a single 
depositional unit. They should also characterize the deposi- 
tional unit, distinguishing it from similar units in a geologic 
basin or province. Within this context, the variables а, 
and ó could represent properties which could be correlated 
in the geologic sense for determining stratigraphic equality. 


If data were available in sufficient quantity to yield 
o 
С. 
Э 


and 6 parameters for individual wells, then these 
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variables could be contoured for specific reservoirs or 
used to segment the reservoir into smaller more homogeneous 
units for mathematical analysis. The values themselves 
would then provide the information necessary for preparation 
of a mathematical model of the reservoir including the 
effects of heterogeneity. 

Thus in reservoir models there is no "a priori" reason 
why the Gaussian normal distribution function need be used 
to study the reservoir performance. If experimental evidence 
indicates the existence of a NT Pearson type 
Orstribution function for a particular property then this 
specific distribution function could be used. 

However before quantitative use is made of the 
statistics а, and ó, it should be emphasized that these 
apply only to the set of data from which they were extracted. 


They are only estimates of the parameters ost 


and ó where 
r r 


the subscript r denotes the value of parameters for the 
population of samples which comprise the entire reservoir. 
For the normal distribution which is described by the mean 
and standard deviation, techniques are well known that 
employ these sample statistics, X and o, to estimate 
population parameters. However very little information is 
available for estimating population parameters а. and о. 


from the sample statistics ET and о. 





CHAPTER IV 
‚ANALYSIS OF LOGRITHMS OF DATA 


Шао 

Jan raw? in his work with core data concluded that | 
"with some exceptions permeability....assemblages give 
satisfactory logrithmic....normal frequency distributions." 
To investigate the applicability of Jan Law's conclusion to 
the core data available for the fields under study is the 
purpose or this chapter. 

It can be observed from Figures Y and 8 that the 
frequency distributions of the permeability of Fields 1 
and 2 are markedly skewed. The distributions are of the 
Pearson Type lj. Does the distribution of the logrithms of 


these data approximate the normal curve? 


Logrithmic Distributions 


If the distributions of these data do satisfy a 
logrithmic normal frequency distribution, an accumulated 
frequency curve of the permeability data plotted on logrithmic 
probability paper should follow a straight line. | 

The sample data may be cumulated and put in percentage 
form as in Tables V and VI. These cumulative percentages 
may then be plotted on logrithmic probability paper. If 
the resulting curve is approximately a straight line, we may 
proceed with assurance to fit a normal curve to the logrithms 
of the data. 


Three cycle logrithmic probability paper was used to 


plot the cumulative percentages from Tables V and VI. 


مم 
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TABLE V Cumulative Distribution of Permeability for Field 1 


Permeability in md. Number of Measurements Percent of Total 


.less than 1.60 120 10 
2.89 195 15 
6.25 260 20 

11.00 325 25 

16.00 _ 390 30 

24.00 155 55 

50.00 520 LO 

37.00 585 45 

45.00 651 50 

53.00 15 55 

65. 00 780 60 

71-00 845 65 

34.00 910 TO 

97.00 5 15 

120.00 1040 80 

126.50 1105 .. 85 

151.50 ITO 90 

185.00 1235 95 

159.00 1505 100 


TABLE VI Cumulative Distribution of Permeability for Field 2 


Permeability in md. Number of Measurements Percent of Total 


less than Jes Da 5 
| 1.8 65 10 
DT .94 15 

[О 126 20 

10.0 158 25 

15.0 189 30 

20.0 220 35 

24,9 252 LO 

29.0 283 45 

52-0 E 50 

36.0 346 55 

2.0 578 60 

48.0 К 109 65 

53.0 447 0 

65.0 172 75 

1220 504 80 

90.0 DOE 05 

122568 572 90 

15.9 605 9 

380.0 635 100 
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The resultant curves obtained for the two fields are 
shown in Figures 24 and 25 respectively. From visual 
observation it appears that the two curves do not approxi- 
mate straight lines. Therefore it may be concluded that 
the permeability data for these fields do not follow 
satisfactorily the logrithmic normal distribution. 

The question that now arises DS LC logrithms of the 
data are not normally distributed, what distribution pattern 
do they follow? To answer this question, the permeability 
data were converted to logrithms. 

The logrithms of the permeability of the two fields 
were then treated as data to which a Pearson frequency 
curve was fitted. The calculated results using the logrithms 
are shown on the following pages, with the calculations and 


figures for Fields 1 and 2 respectively shown. 


Fitting the Logrithmic Data to Pearson Curves 


The procedure outlined in Chapter III was used to fit 
the logrithms of the permeability of Fields 1 and 2 to a 
Pearson type curve. The curves obtained are shown in Eigures 
26 and 27 and were found to be Type 1, and Ig. Additionally 
a normal curve fitted to the histogram is shown for compara- 
tive purposes. 

To simplify the interpretation of Figures 26 and 27, the 
abscissa of the curves were converted to an interval scale 
which is denoted as a m scale. Tables VII and VIII show the 
logrithm value and the millidarcy value of each interval of 


the abscissa. 
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Calculations log K--Field 1 


Data Range (-1.30103, 2.64246) Mean og = 1. 4288151 
Class Interval = .30000 с = „74250 
E = „86963292 & = -, 47453376 Type I, Curve 
гу = -1. 0619523 гр = 5.0271256 
п] = -. 12184190 m, = . 63950490 
С = 7.5821862 
Y = 7. 5821862 (%5+1.0619525)-* "281190 (5,0271256-&)*62950490 
Range] op LEI то oT) 
17 Curve Normal Curve 
Х2 = 50. 551208 X? = 428.38022 
Mic is good at .O0l level Fig 15 пос боса 
т 5са1е Graduation Graduation Normal Curve 
Midpoint Frequency wig-ordinates Areas Mid-ordinates 
1.5 17 сое 226 6.1 
2.5 53 38.5 38.5 16.1 
25 59 poms DTO 16 
1,5 үл 68.6 Og ТЕ) 
5.5 57 85.5 85.6 121.3 
6.5 89 105. 4 105.6 170.5 
7.5 132 130.6 15100 203.8 
8.5 209 16637 OS 206.9 
9.5 297 5T PON TIT S o. 
10.5 266 508.8 289.1 150.9 
17.5 49 or O50 81.6 
12.5 1 0. 0.0 45,2 
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Calculations log K--Field 2 


Data Range, og (-1.2040, 2.5800) Mean og = 1.373706 
Class Intervals = .4000 о = „6920 

а, = .625648 ö = -. 2210993 Type Ip Curve 
ry = 1. 5577193 г. = oo pills 

m, - .62525950 m, = 4.3013559 


8.9991354 x 10”° 


(8.999135!x1072) (t41. 5577193): 09325930. (5. 9815713-«) * 2013559 


У = 
Range] op (-2.14621, 2.473790) 
Lg Curve Normal Curve 
X? = 24.989362 Х2 = 108. 1752 

Fit is good at .005 level Fit is not good 
T Scale Graduation Graduation Normal Curve 
Mid-point Frequency wi g-ordinates Areas (mid-ordinates) 

ШЕ 9 Т 7.4 245 

2:9 20 ET 18.0 10.9 

DES 52 р Som 34.5 

1.5 ЭТ SYe M 626 [io A 

55 76 96.1 96.2 126,3 

6.5 155 DR Тот. = 146.5 

n 178 1559 155.1 Teer 

8.5 90 130.1 ocu (от 

9.5 5 0.0 9.9 ENS 
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Summary 

From the cumulative frequency curves and the logrithmic 
Trequency distribution curves it may be concluded that the 
permeability data under study does not satisfactorily fit a 
Borrichmie normal curve. It is not the intent of this study 
to dispute Jan Law's conclusions that the logrithms of permea- 
bility data approximate a normal curve, but more to express 
a word of caution in indiscrimately using the logrithmic 
normal curve to analyze the permeability distribution of a 
particular set of data. The data used by Jan Law came from 
a different geographical area and from fields with different 
depositional properties than the data for this study. 
Therefore it appears that each and every set of data must 
be analyzed to determine its applicability to the logrithmic 
normal curve. 

Logrithms do provide a means to remove some skewness 
Троп a set of data, but not necessarily to convert it to a 
normal. distribution. If sampling techniques were developed 
for Pearson's Туре lg distributions, the conversion of 
permeability data to logrithms for analysis could serve a 
very useful purpose. At present though, this study indicates 
that caution must be exercised in placing too great of a 
degree of reliability on conclusions drawn from the analysis 


of the distributions of the logrithms of permeability data. 





CHAPTER V 
ANALYSIS OF VARIANCE 


Envrroduction 
From each core analysis there may be readily obtained 


a sample mean (X) and a sample variance ( where: 


(5-1) Kell "a 
and, 
=>, 2 
Z(X, - X) 
=. 
(5-2) Oy = n " 
where X = individual measurements within each well sample 


n = number of samples. 

Within a field the various means obtained from the 
different wells most likely will not be numerically identical. 
It 18 desired to determine whether or not the difference 
between the various means, X, X, 55 x can be explained 
5 random errors, that is, the means are statistically the 
Same and represent identical populations, or whether the 
means are actually different where the differences do not 
ено, Lrom sampling errors. 

To illustrate the techniques involved in making these 
determinations, the data obtained from the core analysis for 
porosity of Field 2 will be used. From this data the follow- 
ing sample means and sample variances were calculated for 


Nude 


the wells cored: 
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Well Number of Mean (X5) Variance (Sy 2) 
Measurements ф 

1 2l . 199 . 000198 

2 28 .189 .000202 

2 go . 200 . 000128 

u 27 .205 .000255 

5 18 .199 . 000263 

6 36 . 186 . 000890 

T 14 . 18H . 000015 

8 25 . 208 . 00017! 

9 34 "ES . 000607 
ТО 20 . 190 . 000228 
mea E 1205 ‚000402 
12 | 19 . 208 . 000894 
13 15 1188 . 000946 
14 _30 қ ‚212 . 000512 

345 2. (15 


From the above it can be observed that the means range 
from .181 to .212. The question which arises, is whether 
ordinary random sampling errors account for the differences 
in these means, or may it be concluded that the means are 


different because of reasons other than sampling fluctuations? 


The F Test 

The F test devised by R. A. Fisher? and named in his 
honor is a means for answering the above question. The 
Basis СОР this test lies in the availability of two independent 
estimates of the population. variance. Consider a single 
core analysis Wló consisting of all of the porosity measure- 
ments from Well 1. The variance Sung of that set is a 


+ 


measure of the internal scatter of the measurements of the 
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population Similar statements apply to S 2 5 Е : 
• \2%? Wo оо о ө Sune 


However, these variances are influenced by sampling error 
md are not necessarily equal to the population variance. 

A better estimate of the internal scatter of the measurements 
may be obtained by pooling the individual-sample variances. 
The better estimate of the within-sample or error variance 


is given by: © 


n n 
L => Ye 2 у \2 "к E 
2 X OA) *X (Хур-Хр) +... 2 (Ху Sine) 


(°) Se 

(n4 + п, Ет 
where na denotes the number of wells and M3» Mo +... Ny 
denotes the number of measurements in wells Wl, W2 .... Wk. 


The numerator of equation (5-3) is called the "within- 
groups sum of the squares", and дес may be called the mean 
Square of individual measurements. 


The degrees of freedom associated with зе“ аге: 2 
т 
(5-4) је = 2 (nk), - п; 


where n, denotes the number of wells and (пк). denotes the 


і 

number of measurements in the ith weil, or fe = the total 

number of measurements in all wells less the number of wells. 
Another independent estimate of the population variance 


may be obtained from the sample means. The variance of the 


population of means is:? ^ - 

5 1 

2 DC je 
(5-5) m, D = n, ==] 
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where: 


(5-6) Mr I + Dot. Ir 


р Number of Wells 
However, the variance of the population of means is not in 
itself equal to the population variance. If all samples 
were of the same size, the second estimate of the population 


2 


variance would be: 
(5-7) S = n S 


If the samples are of different size, n. £ Ny Á thes 


which is usually the case in core analysis, a pooled result 


2 


is to be used as an average sample size, 


n, 

à ol PM 

E aa (2 i 
y 





Therefore: 


(5-9) s ° - no 52 m,p 


? 1s estimated by either equation (5-7) or (5-9), 


2 


wh > 
еп D 


n=l degrees of freedom are involved. is often referred 


to as the mean square of sample means. 
If the random factors which give rise to the within- 
sample or error variance Se“ are the only factors causing 


the differences between the sample means, The Two independent 


2 2 


variance estimates Se^ and Sg should, except for sampling 


error, be equal. The probability of particular rations of 


2 2 


So to Se^ has been computed by Fisher. This distribution 
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is the well known F du M which is a function of the 
degrees of freedom associated with both S ^ and Se*. There- 


fore the ratio, F,’ 2 
5 


(5-10) F = В. 
Se 


is a measure of whether or not random sampling error can 
account for the observed differences between sample means. 
By the use of the F table and the degrees of freedom associated 


E and Sef, a determination of whether random sampling 


Watch 5 
P в 
errors could account for the differences between the means 
can be made. 
conducting the F test on the porosity core analysis 


results for Field 2 gives: 


n n 
1 2 2 = 
2 (X =» 199) + > (X "e 189) + e o e o o 
(24 + 28 о о ө ө ө ) = jJ 
де“ = 
n. 
Y SE 57 2 
2 її 
(5-5) 5 m,p = n. = 


<l 


_ .199 + .189 + .200 .... + .212 
(5-6) DT ee > 





E => 
= 0.128 х 10 
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2 e 
(5-9) 5 
12 + 282 + 268 .... + оё 
(5-8) n = 51-7 и 
о 14-1 224 4-28 + 26 .... + 50 
n, = 24.. 507 
Therefore: 
S ^ = (24.507)(9.128 х 10-2) = .002257 
g 2 
(5-10) pe – „99222! 
| Se . 0004456 
F = 4,907 


Since 5, is associated with 14 - 1 = 13 degrees of 
freedom and Se^ with 345 - 14 = 331 degrees of freedom, the 
Е PUE gives F = 1.76 at the 5 percent confidence level, 

F = 2.14 at the 1 percent confidence level. 

The F value of 2.14 at the 1 percent level indicates 
that there is only 1 chance out of 100 that the observed 
differences between means can be explained by random errors 
if the calculated F exceeds 2.14. In this case 4.907 > 2.14, 
therefore it is reasonable to conclude that the populations 
represented by the sample means are actually different and 
the difference cannot be explained by random variation. 

Since there is an actual difference between the 14 
Sample means, further testing must now be conducted to 
attempt to determine if there are sub-groupings of these 


means that are statistically homogeneous. For this field, 
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the wells were lumped together by leases and an F test 
conducted to determine if the means of the three leases 


were statistically the same or whether they were different. 


Field 2 
Lease Number of Mean (X¢) Variance (5,2) 
Measurements 
all 179 .1948 . 000383 
va .1950 .000670 
3 _95 20р ‚000641 
345 
For leases: 
бес = 000512 
2 
Sn = „002368 


с 2 
F = 5 = 1.62 
se 


Degrees of freedom with зе“ = 545 - 5 = 342. 
Degrees of freedom with 5,4 = m 
The F table gives: 


3.03 at the 5 percent confidence level 
4,69 at the 1 percent confidence level 


Since the calculated F, 4.62 < 4.69, the observed 
difference in the three leases means could reasonably be 
explained on the basis of the scatter of the observed data, 
and the three means are statistically the same. 

It may be concluded from the two F tests that the 
individual well means are statistically different, but the 


lease means are statistically the same. Ihe well means each 
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represent individual points in the field and since these 
points are different, heterogeneity exists between the 
various wells. As the wells are lumped together to form a 
lease mean the heterogeneity of the individual points are 
absorbed into larger somewhat homogeneous units, these 


larger units exhibiting similar statistical characteristics. 


Modified Tukey Test 

The results of the F test on the means of the 14 wells 
gave convincing evidence of differences among the means, 
but the F test gave no clue as to how many differences there 
were. Ina group of a means there are in all a(a-1)/2 
potential differences; 14(13)/2 = 91 among the wells. Does 
each mean differ from all the rest, or are some of them 
the same? 

One method of investigating the differences is by the 
Tukey test (modified)?. The test is made by computing a 
difference, D, which is significant at the 5 percent level, 
then comparing it with the a(a-1)/2 sample differences. D 
Ше рас product of ox and a factor, Q, taken from a Q table? 


which is itself computed on the basis of the Оле Био ош 


the deviations among compared means. 


© 
X 
5-11 Ss = == 
(5-11) к 
and: 
(5-12) D= Sy 





ШІ 
То determine Q enter the upper heading of the Q table 
with the number of treatments (14 wells) and degrees of 
freedom, f, for samples (331) indicated at the left of the 
table. From the table with 14 treatments and f = 331, Q 


is 4.74. 
BEES oon, Field 2, 

Sz = 1.86 x 107° 
and, 

Se = 4.315 x 107 
Therefore: 

D = (4.74)(4.313 x 1077) = .02048 
The differences to be compared with D are shown in 

Table IX. 


In Table IX the X, are arranged from high to low and each 
is subtracted from those above. Of the 91 differences, only 
15, indicated by *, exceed D = .0205. One inference from 
the leise that the population represented by the mean of 
well 13 is different Pon that of wells 14, 8, 12, 4 and 11. 
Similar statements may be made about any particular well. For 
instance it may be inferred that the population represented 
by the mean of well 14 is different from that of wells 13, 7, 
6, 2, 10, and 9 but the differences between well 14 and the 
remaining wells are not significant. This above procedure is 
therefore useful for seeking homogeneities and dissimilarities 
among -wells. Such information could be valuable in explaining 
Similarities and differences in individual well production 


behavior. 
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sequential Method of Testing 
A sequential method of testing the differences between 


7 


the means devised by Hartley is a somewhat more powerful 
testing procedure than the modified Tukey method. For this 
test, not one Q is taken from the Q table but several, one 
for each range of the treatment means. For the well means, 
adjacent means in the array are tested with Q = 2.77 for a = 
2; for two ranks apart in the array use 9 = 5.52 Гог а = э; 
for three ranks apart in the array use Q = 3.63 fora = 4; 


use Q = 4.74 only for the extreme range where a = 14. The 


corresponding D are: 


а, 9 D 
2 А ROT 
3 nuo kom 
П 3.63 . 6 
5 3,86 Bo: 
6 1.03 O14 
7 LT O16 
8 1.29 . 0185 
9 1.39 . 0189 
10 4.47 , 0193 
nen 1. 55 · 0196 
12 102 .0201 
13 к бо оор 
14 ye . 0205 


These D apre entered in the northeast-southwest diagonals 
of the table of differences, Table X; with the D's in 
parenthesis. Each difference is now compared T ICS Ош Ш 
the difference being judged significant if it is larger than 
its D. A useful rule to be observed is that if any difference 
is less than its D then no further testing needs to be done 
to the right of that difference in its row or below it in 


CE COLUMN. 
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From Table X, it is observed that there are 18 signi- 
ficant differences, marked by an *, which signifies that the 
sequential method detects a greater number of differences 
than the modified Tukey method. 

From the table the wells similar to any given well may 
be detected. For example, it may be stated that the mean of 
well 9 is not significantly different from any of the other 
means except well 14. The mean porosity of well 3 is not 
Significantly different from that of any other mean except 
for well 13. 

The sequential testing method then provides a comparison 
mechanism for determining the relationship between any well 


mean desired with the other well means in the group. 


Homogeneity of Variances 


Several tests are available for testing the homogeneity 
of the variances of the several samples. The appropriate 
test depends upon the type of possible variation of the 


variances that is visualized. If the deviation from equality, 


ea (a, = T I ME o, *) is conceived to be caused by 


a random variation then an appropriate variance homogeneity 
cest is Bartlett's — wherein the statistics B and C are 


computed where: 


1 252 | 2 
В = у(у Ја 5 — 2 di Se) 


(5-13) 


а = 


EDO 209 2 y 
= + (v 1051 0 S Z у, logig 81^) 





LTO 


s _i 
Va V 
(5-14) C=1+ 
5 (K-1) 
wnere: 5 
SVS. 
Vv = 2 v. and 52 E E 
i V 


Va is the degrees of freedom associated with each sample 
variance БЕ and K is the number of variances being con- 
Sidered. 

The statistic B is known to satisfactorily approximate 
the Chi-square (x=) distribution corresponding to K-1 degrees 
of freedom.” A calculated B value greater than the x° value 
at the given degrees of freedom and confidence coefficient 
is evidence for rejecting the hypothesis A zc Ga dC us 


: E. Calculations: for porosities of Field 2 are given below 


to illustrate the method. 


TABLE XI 


Computation of Bartlett's Test of 
Homogeneity of Variance-Porosity Field 2 


Well "Vi 2 * 2 2 2 2 

No. (n-1) 3i m E m T Vi Уі 94 
al 23 ко зоо 6.8384 ‚044873 1.5.609 
2 оү 0200 OS 1 9.2599 057057 54,513 
3 25 о 2,7058 . 040000 2.078 
Д В 551 29967 10.3914 ‚03846 66.526 
5 17 о №57 ЭЛЭ . 05882 т 
6 35 8.903 „94954 29. IG Eoo 511.605 
7 15 ir 02.111950 TING OOD 106,79 
8 22 1.720. 24055 5. 2921 . 94546 35.250 
9 33 SUOMI 76520 25.8476 . 03030 200.545 
10 19 UE SO) M 6.8081 · 05265 15.555 
d 30 1.025 „60155 18.1365 205225 1202098 
T2 18 8.938 .9512! 17.2225 . 05555 160.884 
13 14 9.458 .97530 13.6612 . 07145 152.412 
11 29 35.109 .H9262 14. 2860 . 03448 OON 
331 TUAE EO .64788 1557.815 


*Converting porosity data to percent (.199 to 19.9) changes 
variance from .0001985 to 1.985 percent. 





(5-12) i 
(5-1!) С 
| 
C 
x 
C 
(5-12) I 


UJ 
|| 


UJ 
|| 


From ж table” 


K-1F v-221 A 


а-20224 (v los S“ - >x V4 108 Ta 
т 

~ 

A о 
5 (x1) У 2M 
_ 1 - „84788 _- „0030211 
3 (14-1) 
=1+ .01153 = 1.01433 
2 

AS E 

eE M i = 4,1022 


1987 0 S^ = ‚61302 


(02252) ( (551) (.61202) - 171.2596) 
(2.27006) (51.67) 


1:99 


with 13 degrees of freedom (14-1) 


а = 22,4 (at 5 percent level) 
хе = 27.7 (at 1 percent level) 
E = 29.8 (at 1/2 percent level) 


A у“ of 29.8 at the 1/2 percent level indicates that 


chere is only 


1 chance out of 200 that the differences 


between the variances could be caused by random variation if 
M. 


the calculated B » 29.8. Since the B = 71.98, it may be 


coneluded that the differences between the variances of the 


14+ wells are significant and cannot be accounted for by 
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random variations. Therefore the hypothesis that ЈЕ = 
2 oe 

р =.... O0. 15 rejected. 

summary 


Various statistical tests that may be used to analyze 
core samples by use of the core means and variances have 
been illustrated. A core sample provides certain information 
about the population from which it is extracted. The res- 


pective populations sampled by core,, core ... соге- у, 


2 

fOr Fiela 2 may be thought of as falling into one of four 

categories: 

1. The 14 samples represent identical populations, the 
apparent differences resulting from sampling error. 

2. They represent identical variances but have different 


means. 


They exhibit identical means but different variances. 


Vl 


They are totally different. 


The F test indicated that the populations had different 
means. Bartlett's test showed that the populations had 
different variances. Therefore we may conclude that the 
populations represented by the 14 samples are different. 

By use of the Tukey test, and the Sequential test it 
is possible to determine what ‘sub-groupings of the 14 
populations are statistically homogeneous. An indication 
that any particular well is statistically the same as a 
sroup of other wells should prove to be a valuable tool in 


correlating the expected performance of that well with the 
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performance of the group with which it may be identified. 
| In this chapter particular attention has been given to 
the distribution of sample means and of sample variances, 
yet in reality there is very little exact knowledge, based 
upon sedimentation theory, as to the appropriate distributions 
to expect. It is to be hoped that the methods of the 
present chapter will be useful for quantitatively character- 
izing these distributions and thereby aid in the formulation 


of more precise theories of sedimentation. 





CHAPTER VI 
SAMPLING 


EurpPoduction 

The prime purpose of core analysis is to gain some 
understanding of the characteristics of the field properties 
under consideration, namely porosity, permeability, and 
Fluid saturations. Ав previously discussed, once sufficient 
information becomes available, it may be used to determine 
descriptive parameters for a field. The parameters may be 
expressed simply in terms of averages and deviations from 
the averages. Or more extensive parameters may be developed 
ШІ ters Ol theoretical frequency curves anid the moment 
parameters апа д. 

wereracring the initial period in The ITE GOT Е еле 
before an extensive amount of data becomes available, the 
question arises as to what inferences can be derived from 


the data as it becomes available. 


sampling 

Sampling may be defined as the selection of part of an 
aggregate or population, on the basis of which a judgement 
or inference about the aggregate or population is made. Y 
This sampling theory is a study of relationships existing 
between a population and samples drawn from the population. 
Statisticians have long worked within the problem of reconstruc- 
tion of a universe of variables by means of samples that 


comprise a small percentage of the universal or population 


from which the samples were drawn. Core analysis data appears 
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to comply with the requisite of random sampling as stipulated 
ES Theoretical statistics. 

| Sampling theory is useful in estimation of the unknown 
population quantities such as population mean, standard 
deviation, variance, etc., referred to as population para- 
meters, from a knowledge of corresponding sample quantities. 

In general a study of inferences made concerning a 
population by use of samples drawn from it together with 
indications of the accuracy of such inferences using prob- 
ability theory, is called statistical inference. 

Considering a core analysis as a random sample drawn 
from an infinite population, the sub-surface of the earth, 
what inferences can be made about the population? From the 
core analysis there may be obtained a sample mean and a 


sample standard deviation where: 


= XX 
(6-1) X = = 
И 
(6-2) کلک = رو‎ 


The sample mean, which serves as an estimate of the 
population mean will be of chief interest in considering a 
core sample. An important theorem in statistics. states 
пас for almost all populations the probability distribution 
or the sample mean based upon a simple random sample will be 
an approximately normal one if the sample size is sufficiently 
large. "© The standard deviation of the probability distribu- 


tion of the sample mean, denoted by Or is calculated ъу:2 
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(6-3) OF = 


5 | 


Where Oy is the standard deviation of the population being 
sampled, and n is the sample size. 

Therefore if the standard deviation of the population 
is known, Oy nay be calculated. However, after only one 
core analysis, the standard deviation of the population is 
not known and an estimate of this standard deviation must be 
made. А point estimate of this population characteristic 
is made from the sample standard deviation, Sys and the 
Be mean is the best point estimate available for the 
population mean. 

To illustrate the application of point estimates, the 
numerical data obtained from an assumed first core analysis 
taken from Field 2, Figure 28, will be used. 


From the first core analysis, Well 1, we obtain: 


Ky = 18.6% 
Syg = 2.98 
Хо = 15.82 
5 от №56 
Х - 15.55 
Sy 14 

n = 36 


Therefore the best estimate of the various population 


means and standard deviations are: 
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FIGURE 28 MAP OF WELLS CORED--FIELD 2 
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=> Е c 
^ópopulation С 127 
C ` = 45,8% 


SOpopulation 


- 45.37 


SW population : 


ра] 


and, 


And, we assume that: 


Oz = — — = : 
s Ма 3 

go = ‚70 

Fay = +69 


For normal probability distributions, the area under 
the normal distribution curve between the mede Standard 
Y 


deviations is about .95 out of a total area of 1. Therefore 


the following limits may be constructed with a confidence 


on. 95. 

(6-4) X+20X 
X, + 2(.48) = 18.6 + .96 

17.64 - 19.56 

2(.76) = 45.8 


14,2 - 47,4 


Ра 
8 
1+ 


sot 1.52 


Ра 
er 


2(.69) = 45.3 + 1.38 
43.9 - 46.7 


ә = 
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We can conclude with a .95' probability that: 
^o 
o for Field 2 is somewhere between 44.2 and 47.4%. 


for Field 2 is somewhere between 17.64 and 19.56%. 
E. for Field 2 is somewhere between 13,9 and 46.7%. 

The above conclusions are based on knowledge of the 
true value or oy, which was not the case. Statistical 
history has indicated though that if the sample size is 
reasonably large, n > 30, the Sy may be used as an estimate 
о? Oy without materially changing the reliability of the 


above conclusions. 


moncrol Charts 


From the previous section, estimates of the probable 
range of the various property means of Field 2 were deter- 
mined. The ranges were based upon one cope analysis which 
Was taken from one point in the field. For the population 
which this core sample represents, the sample mean ranges 
determined should prove to be satisfactory range estimates 
of the population mean. As drilling moves away from this 
point, a different areal population may be encountered with 
different characteristics from those determined by the 
initial core sample. In this case it would be expected 
that the new population characteristics may deviate from 
those estimated from the initial core sample. 

In statistical quality control procedures for production 
processes, a control chart is used to indicate when a process 


has changed or "gone out of control" and is no longer 
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producing within prescribed specifications. Basically the 
production control procedures consist of computing a process 
mean and standard deviation, constructing a control chart 
within the limits of Xprocess + 3 Oy, taking periodic 
samples from the production line and computing the mean of 
the samples taken. If the mean falls between the limits 
set on tne control chart, it is concluded that the process 
has not changed or “is in control." If the mean does not 
fall within the limits, it is concluded that the process has 
changed or "is out of control. "2 

Core sampling may be considered analogous to production 
line sampling with areal sampling of populations considered 
tO be the analog or periodic process ER As core 
samples are taken at varying distances from the initial 
sample, the means may be plotted on a control chart to indicate 
if the population: characteristics are still "under control," 
or if not, whether a new areal population has been encountered. 
For a control chart, + 3 o= is normally used. Thus: (from 


A 
previous section) 


5 обу = 18.6 + 3(.48) 


X, £ = 
18.6 + 1.44 
17.16 - 20.04 
212 сло M 
45.8 + 2.28 
43.52 - 48.08 
x 2 05 = 45.3 № 21-69) 


SW 


. [T 
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Using the above, a control chart may be constructed. 
From Figure 28 it can be observed that the wells are located 
in a general north-south direction with well 1 located 
approximately at the center of the field. 

lun this particular geographical arrangement, two 
control charts to test the hypothesis that the mean of the 
field falls within the probable range are constructed, one 
for moving in a northerly direction, one for the southerly 


direction. 


For moving North: 


To 


2.1 


(Upper “A 
Limit 20.04 


20 


Average 18.6 


ë (Lower Control ) 
о ЕКЕ 


17 


NI A2 A3 NY NE N6 
sx. Wells 


FIGURE 29 Control Chart for Average Porosity 
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The core analysis from well Nl has а Xo 


is plotted on the control chart, and falls within the limits 


of 19.2. ers 


set. It is concluded tnat the sampling is still under 

control T tne population has not changed from that indicated 
by the initial sample. As the core analysis from well N2 
through NÓ are taken, their means are plotted on Figure 29. 


Moving north from initial well: 


Well Xo 

№ 19.2 under control 

N2 20.5 POSSIDLY outo Coni oi 
№5 19.9 under control 

N} 10,9 under control 

N5 19.9 under control 

NÓ 19 O under control 


It appears that the probable range of the field mean 
calculated from the initial core analysis gives a good 
representation of the field mean from the initial point to 


well N6. 


For moving South: 


о 
21 
(Upper control) 
Limit 20.04 


ёо 


I3 
Average 18.6 


/8 
(Lower Control, 
a ШЕ 


/6 
Si $2 $3 59 56 


Wells 


FIGURE 30 Control Chart for Average Porosity 
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Moving south from initial well: 





Well Хе 
E qut under control 
52 19.1 under control 
S3 20.8 possibly out of control 
SA 20.8 more indication of 
OUT Of сопсто 
95 21.2 out of control 


For movement in a southerly direction the initial 
conclusion concerning the field mean proved correct from the 
initial point to well 532. When it becomes apparent that 
the initial hypothesis concerning the mean is not tenable, 
in this case after well S4 was drilled, a new probable range 
of the areal mean should be calculated from the data obtained 
from well S3 and S4. This new probable range would then be 
the best estimate of the population mean for the areal 
population from S4 southward. 

The adoption of the above method, which is based on 
production quality control techniques, to the analysis of 
reservoirs is one possible way in which a quantitative criterion 
could be used to segment a field into smaller homogeneous 
units. Many extensions to this procedure are possible. 

These involve the comparison of the segmentation created by 
use of the different reservoir properties. In this regard 

it is possible that the areal segments would not coincide 

when using different properties such as porosity, permeability, 
or fluid saturations. However there is also no "a priori" 
reason for their being different. Similarity of areas 

would probably indicate a fairly high degree of correlation 


between the variables involved. 
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Confidence Interval for Population Proportion 

If a random sample size is large, n > 30 usually being 
considered sufficiently large, a confidence interval for the 
population proportion can be constructed in a similar way 
as for the population mean, since the probability distribu- 
tion of the Sample proportion will be assumed to be approxi- 
mately normal. 


ithe standard deviation of the probability distribution 


of the sample proportion can be estimated on 
(6-5) | L 
è n- 


where p is the sample proportion and 55 denotes the estimate 
of the true standard deviation ex 


The confidence interval for the population proportion 


where Z is the normal deviate corresponding to the desired 
confidence coefficient. 

To illustrate the М” of confidence intervals 
for population proportions to core analysis, suppose that 
it has been decided that an oil saturation of 50 percent or 


less is undesirable in considering the amount of pay zone in 


a field. The initial core analysis taken from Field 2 will 


be used for example calculations. 


There were 36 samples from well 1, 6 of which were < 30 


percense,.. Thus: 
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and, 


s= = Dee ed ЕЕ сс 


And tor a confidence coefficient of .95, Z is 1.96, ° 


u 
EID 


| 


.167 + 1.96 (.066) 
.167 + .129 


Thus it may be concluded with a confidence coefficient 
of 95 percent that the percentage of the field pay zone that 
is undesirable according to the predetermined criteria is 
somewhere between 3.8 percent and 29.6 percent. 

Another way to interpret the use of the population 
proportion and its relation to the confidence coefficient is 
illustrated by the following modification of the above 
problem. Instead of establishing a confidence coefficient 
oi 0.95 and determining the range in the percent of formation 
which would be unproductive, suppose that one Seeks to know 
the confidence coefficient associated with a particular 
percentage range such as 1/6 + 3 percent, or .167 + .03. 

One then determines the value of Z corresponding to cne T 


2 percent range. In this case: 


interval in population proportion =p + Z S= 


A value of Z is sought such that 2(5-) - .05 ог 2.- 
.03/.066 = . 455. This value of Z corresponds to a confidence 
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evefricient of .35. The interpretation is that the true 
unproductive fraction of the formation has 2 chances out of 
5 to be outside of the proportion .167 t TOJ or a menange 


in three (approximately 35 percent) of lying within the 


proportion range .167 + .03. 


Summary 


Two rather elementary techniques of estimating informa- 
tion about a field waen the amount of core analysis data is 
limited, in this case one well, have been presented. For 
Field 2, the sample mean of well 1 gave good results as an 
estimation of the population mean for the areal population 
from well S2 to NO. 

Tne utility of these sampling techniques can only be 
confirmed after being tried and tested with a large number 
of fields where the actual results could be compared with 
the estimated predictions. It would be conceivable to assume 
that over the long run, sampling techniques would provide a 
better estimate of what is to be expected than those estimates 
determined intuitively. Granted, in some cases the reservoir 
predictions by sampling techniques would prove to be in error, 
Се иоле ево штасев could show less error then vhezerror 
or difference found by intuitive predictions. 

Sampling results SEM be recognized as only estimates 
and not as infalliable predictions. Yet the sampling 
estimates should prove better than no estimates at all. 


For any one reservoir the estimates are either right or 
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wrong. The .95 confidence coefficient means that if the 
procedure is followed for a long enough time over a sufficiently 
large number of reservoirs, the results should be correct 
95 percent of the time. 

It was the purpose of this discussion on sampling To 
indicate how certain techniques could be applied and possibly 
poostumulate the study of the applicability of other more 
advanced or specific sampling techniques. It is felt that 
the work and investigation in this area has just begun with 


unlimited possibilities yet to be explored. 





CHAPTER VII 
SUMMARY AND CONCLUSIONS 


Introduction 


At the outset of this study two alternate courses of 
investigation none considered ae the statistical analysis 
of the field data available; the first, to concentrate on 
one Statistical technique and analyze it in detail, and 
Second, to explore various techniques with less detailed 
analysis oi any one point. The latter course was chosen 
with three very broad areas included in the study, i.e., 
generalized or skewed distribution curves including trans- 
formation of data to logrithms for analysis, analysis of 
variances, and sampling procedures. 

The choice of this latter course was prompted by the 
realization that relatively little use of theoretical 
Statistics has been made in core analysis as judged by 
available discussion in the petroleum literature. An 
attempt was made to present in one volume the necessary 
calculation procedures and interpretations to serve as a 
broad guide for & statistical analysis of field data. 
Additionally, it is hoped that the study will serve as a 
basis for further investigation into the applicability of 


the science of statistics to reservoir problems. 


Summary 
This thesis has attempted to indicate the applicability 


of certain statistical procedures to the analysis of core 


samples. Statistical techniques were employed to segregate 
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relevant and meaningful information from a mass of raw data. 

A method of determining and describing frequency distri- 
butions of the physical properties of a reservoir is presented 
in Chapter III. The Pearson family of generalized frequency 
curves provides an effective means of classifying observed 
frequency distributions without being restricted to a few 
classical distributions such as the normal curve, the Poisson 
Emé Or the binomial distribution. 

Chapter IV is devoted to the study of the effects of 
the transformation of raw data to logrithms for analysis. 

The results obtained were compared with the results of Jan 
Law's* work. 

Analysis of variances is Sed in Chapter V with 
various testing procedures shown. The F test and Bartlett's 
test ape employed to test for uniformity of the population 
represented by the various well samples. The modified 
Tukey and Sequential tests are used to determine what sub- 
grouping of the wells are statistically homogeneous. ' 

Sampling procedures are presented in Chapter VI. 

Sample means are analyzed and presented on control charts 
by techniques analogous to production line sampling techni- 
ques. Population proportions are presented with confidence 
intervals. | _ 

The following listing presents a summary of the main 


eguations associated with the above statistical procedures. 





(2-21) 


Type IV Curve 
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For Frequency Distributions 


Pearson Type I Curve 
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poe VI Curve 
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Normal Curve 
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The F Test 
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where: 
(5-9) 


and, 


(5-3) Se 


2 __ 
> ~ (X-X) 
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For Analysis of Variance 
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(n4 +.... n.) zen 


Modified Tukey and Sequential Testing 


I) 
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М... S 
S> = — 
Ja 
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where: 


Et) 


l 
B= (У ln S 


151 
Bartlett's Test 


e . X y, ln 54°) 


1 
:d 

C = 1 + i 

5 (K-1 


For Sampling 
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"Onclusions 


In conjunction with the comments given in the summary 


sections of the specific chapters, the major conclusions to 


be derived from the data studied are believed to be: 


|. 


„ч 


Generalized frequency curves are useful in reflecting the 
variation of individual properties throughout a single 
depositional unit, and in characterizing and distinguish- 
ing the unit from other units in the same geological 
basin. The terms o and ó should provide additional 
relevant numerical parameters for mathematical model 
analysis or a reservoir. 

Logrithmic transformation is occasionally useful to 
remove some skewness from permeability data, but does 

not necessarily convert the data to the normal distribu- 
pop under all circumstances. 

Analysis of variances tests provide a useful tool to 
investigate homogeneity between well properties within a 
reservoir or within sub areas ог а reservoir. 

Sampling techniques may be utilized to estimate unknown 
population quantities from a knowledge of relatively 
small sample quantities. 

The science of theoretical statistics applied to reservoir 
analysis is in its infancy and appears to offer many 


opportunities. Tor the solution of reservoir problems. 
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Recommendations for Turtiier Investigation 


While working with the Pearson frequency curves, certain 


areas of further investigation appeared to offer promise: 


il. 


An investigation into the P between reservoir 
performance and the type of frequency distribution 
iollowed by its physical properties. For example, does 

a field with a Pearson Type III Sy, distribution perform 
differently than a field with a Type 1, other factors 
being the same? That is the correlation, if any, 

between the parameters (с, 6) and field performance, 

the parameters c and ó being considered for any physical 
property or group of properties? 

sampling techniques are available for populations that 
approximate the normal frequency distribution. Could 
reliable techniques be developed for populations that 
satisfy certain types of the Pearson family of frequency 
distributions? A table of Areas Under the Type Ill, 
Curve has been developed” and could possibly serve as a 
point of departure for investigating sample techniques 


for core samples. 


Numerous tests have been devised for the analysis of 


variances. Certain of these tests were presented in this 


thesis. Other tests may prove more useful than the ones 


presented. Further investigation into this area of statisti- 


cal analysis is highly recommended. 
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(са Op., Cit: ED 


Ideen. 21. 
ВО ое 22; 


Craig, OD. Cit., pp. 277247 
Craig in developing the Pearson equations in terms 
of а, and ó, shows this equation but does not 


offer a proof. For purposes of this study, the 
equations for the various constants, C, were used 
as presented by Craig. 


Pearson, Karl, Tables for Statisticians and 


Biometricians, Cambridge: 1914, pp. 120-172. 


ocarborough, JW Numerical Mathematical Analysis, 
Baltimore: 1950, р. 145. 


and 


Lowan, A. N., Bulletin of the American Mathematical 
Society, Oct., 1942, Vol. 35, No. 10, pp. 759-715. 
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Scarboroush shows the method of numerical, integra- 
tion involvins Gaussian coefficients which is 
considered to be the best general purpose quadrature 
procedure. The method for determining I where: 


b 
T G NE 
a 
by this technique is: 


іст 
Т = (Ъ-а) "NS R. £(Z, ) 
wnere m are the selected Gaussian coefficients 


and the Z, are Selected, known values of the 


independent variable x transformed to a scale 
wherein the range a to b becomes the range from 
Otol. For the present integration it was found 
necessary to use n - 16 for appropriate accuracy 
of the function. The "Weight Coefficients" 

for the Gauss Quadrature Formula for order 16 
were taken from Lowan, pp. 739-743. 


mitn, J. G., Elementary Statistics, New York, 
1934, p. 287. 


Carver, OD. отс рош 


ска, ор. сто рс 


Chapter III 


i 


Nl 


When data are grouped in discrete intervals, the 
moments are obtained by a summation process. The 
frequencies are considered to be lumped at the 

mid intervals which may not be the actual case. 
Therefore grouping errors may be introduced. An 
extensive discussion of a method of correcting for 
these grouping errors, known as Sheppard's cor- 
rections, is presented in Kendall, M.G, The 


Пева Theory Of Statistics London, Vol. I, 
1945, p. 71. Sheppard's corrections are: 


X, = Ho - 1/12, X, = Mz, Xy = y ~ 1/2 u, + 7/240 


Kendall, M. G., The Advanced Theory of Statistics, 
London, 1945, Vol. 1, p. 49. 


Ibid, p. 50. 


Davis, H. T., Tables of Higher Mathematical 
Functions, Vol. I, Bloomington, 1933, p. 179. 





C 


10. 


Tl. 


2ا 
Ibid.‏ 
Ibid; pp. 217-239.‏ 


Ibid. p. 181ff, gives this formula and other 
asymtatic expansions of the Gamma Function. 
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Law, Jan, Statistical Approach to the Interstitial 
Heterogeneity of Sand Reservoirs, AIME Transactions, 
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Chapter IV 


ШЕ 


Chapter V 
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Шали Jan, Op. Cit., P: 2215 
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Statistics, Englewood Cliffs, 1959, p. 5. 
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APPENDIX A 


|, MATHEMATICAL, 'FORTRAN, AND TECHNICAL NOTATIONS AND TERMS 


Mathematical Fortran 


Notation Notation 
A A 
ps ASQ 
a 

ДЈ 

АК 

ао 
n 

а ALPHA 
а ALP5 
'ALP3S 
AN 

ARG 

ARG1 


moment above mean 
2 


Description 


A parameter of the Pearson Type III 
Curve, equal to 2/ALP3 


Square of A 


Parameter of the differential equa- 
tion (2-1); also skewness of the 
distribution cf a set of data defined 
by equation (2-15); also, a term used 
in the Sequential method of testing 
to describe the range of means in an 
array 


Number of measurements in a particular 
class interval 


Total number of measurements 


ch 
Then а term where а, = 


„eh 
see equation 


C 


A parameter of the Type VI Pearson 
Curve, see equation (2-37) 


A basic parameter of the Pearson 
system of Frequency Curves, see 
equation (2-7A) 

Square of ALP; 

A pasic parameter of the Pearson 
System of Frequency Curves, see 
equation (2-7A) 

Number of class intervals 


Intermediate variable for computing 
YBAR of the type curve 


Intermediate variable for computing 
YBAR of the type curve 





Mathematical 
Notation 


B 


Fortran 


Notation 


BELSQ 


B5Q 


CHI 


CHISQ 


CON 


CST 


142 


Deseriprion 


A statistic calculatedkin Баро ес = 
test for variance homogeneity, See 
equation (5-13) 


An intermediate value for derivation 
of the constant 2 Type I curve, 
see equation (2-28) 


бо» Ру» b, are terms of equation (2-2), 
where f(t) is expanded in a converg- 
ing power Series. values oí terms are 
defined in equation (2-11 


Summed BSQ between observed data and 
Normal Curve 


2 
Equivalent to an individual 
on page 45 for Normal Curve 


A constant in the frequency curve for 
a particular set of data, also a 
statistic calculated in Bartlett's 
test, see equation (5-14 


= 
Equivalent to an individual . 
on page 75 for Pearson Curve 


Summed CHI between observed data and 
Pearson Curve--a statistical distri- 
bution used for certain statistical 
testing 


Storage locations for all parameters 
and the constant characteristic of a 
given Pearson Type Curve 


Pearsoni O ин m iV 
CON(1 ЕМІ А R 
2 EM2 C EM 
2 RE S 
i R2 V 
5 C C 


The Gaussian coefficients (16) Гог 
integration in determining F(R,V), 
see Ref. (7 Chapter 3 for source of 
coefficients 


An intermediate parameter for all main 
Pearson Type Curves, see equation (2-15) 





Mathematical 


Notation 


D 


>l 


3 


Fortran 


Notation 


DAVG 


DD(I) 
DECM 


DELD 


DELTA 


DMAX 


DMIN 


EM2 
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Description 


For analysis of variance, a computed 
difference between means that is 
significant at a chosen level, See 
equation (5-12) 


Arithmetic mean or average of data, 
see equation (6-11) 


“Midpoint value for class interval I 


Difference between the highest value 
permissable in a given interval and 

the lowest permissable value in the 

next higher interval 


Difference between maximum and 
minimum value of a class interval 


A basic parameter of the Pearson 
System of Frequency Curves, See 
equation (2-12) 


Maximum value for data in highest 
class interval 


Minimum value for data in lowest 
class interval 


Term used to locate interval which 
contains the mean 


The square root of D, see equation (2-15) 


A parameter in the Type IV Curve, See 
equation (2-30) 


A parameter in Type I and VI Curves, 
see equation (2-23) 


A parameter in Type I and VI Curves, 
see equation (2-23) 


The independent normalized variable 
for the interval - 1/2 to + 1/2 based 
upon 16 subdivisions for Gaussian 
integration. The 16 values of EU are 
taken from Ref. |4 Chapter 2 





Mathematical 


Notation 


fe 


Fortran 


Notation 


EZ 


FRV 


G(J) 


ID 


IDEPT 
INDEX 


TE 


IPERM 

IPHI 

ISAMP 
150 


144 
Description 


The term (Sin 0)“ e"? used to determine 
F(R,V), see equation (2-33) 


A ratio devised by Fisher to test for 
homogeneity of means, see equation (5-10) 


Degrees of freedom associated with зе“, 


see equation (5-4) 


F(R, V) function, used for constant 
of Type IV curve, see Ref. ії , Chapter 2 


A function used in the derivation of 
the constant from Type IV curve, see 
equation (2-33) 

Frequency of values in an interval 
Index, usually denoting interval number 
Term to indicate a complex root in 
derivation of Type IV curve, see 
equation (2-30) 


Identification number for the field, 
well, and lease 


Depth of sample, feet 

An index to distinguish permeability 
data cards from other data cards, 
INDEX = 1 for permeability 

INDEX = 2 for other data 

An index to distinguish the first 
cycle of a loop from all subsequent 
cycles 

Permeability, millidarcys 

Porosity, fractional 

sample number 

Oil saturation, fractional 

ТОСА 


lndex denoting number of different 
Pearson Type Curve fits desired 





Mathematical 
Notation 


k 


rol 


Foptran 


Notation 


MP 


NA 


NDEX 


NTYPE 


NUMBR 


Pl 


PROP 
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Description 


Nunber of variances considered in 
Bartlett's test 


Identification Program--numerical 
equivalent of AN in fixed point notation 


Fitting Program--index used in inte- 


‚gration for YAVG 


Number of the interval counting from 
the lowest which contains the mean 


Number of intervals in fixed point 
notation 


An index to distinguish among porosity, 
oil saturation and water saturation 
data cards: ` 


NDE SLU LGreperosi uy 
=)? TOR OIl вариа топ 


5 Рог water saturation 


An estimate of an average sample size, 
see equation (5-8) 


Type of Pearson Curve to be calculated: 


NTYPE Pearson Type 
1 [ 
2 WEF 
2 ТУ 
hc VI 


Number of different type of curves to 
be calculated in a computer run 


An intermediate variable for the 
Type IV curve 


sample proportion, see equations 


A 
(6-5) and (6-6) 


Variable for an individual physical 
property such as porosity, permeability, 
etc. 


An interval. scale for the logrithms 
of permeability data 


A factor obtained from a table to 
compute D, see equation (5-12) 





Mathematical 
Notation 


q 
Aq 


Fortran 


Notation 


R 


R1 


R2 


S 


SIGMA 
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Description 


Mean value of q, see equation (3-5) 


The interval of q used to construct 
the histogram on the q scale 


A parameter of the Pearson Type IV 
Curve, see equation (2-30) 


A parameter of the Pearson Type I and 
VI Curves, see equation (2-15 


A parameter of the Pearson Type I and 
VI Curves, see equation (2-15) 


“A parameter of the Pearson Type IV 


Curve, see equation (2-30) 


The standard deviation of the set of 
data, see equation (3-9) 


A term in Bartlett's test equal to 


DU V а = 
Du ү 


V 


Witnin sample or error variance, 
see equation (5-3) 


sample variance in Bartlett's test, 
see equations (5-15) and (5-1!) 


Variance of the population of means, 
see equation (5-5) 


An estimate of the population variance, 
see equation (5-7) and (5-9) 


Standard deviation of the probability 
distribution of the sample proportion, 
see equation (6-5) 


The sample standard deviation, see 
equation (6-2) 


A sample variance, see equation (5-2) 
The standard deviation of the probability 


distribution of the sample means, see 
equation (5-11) 
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Notation Notation Description 
© The standard deviation of the prob- 


ability distribution of the sample 
mean, see equation (6-3 


SUM An intermediate storage cell used for 
the Gaussian integration in determin- 
ing F(R, V) 


SUMD oum of data values in an interval 


SUMXG(I) Sum for the I'th moment about the 
midpoint of the interval containing 
the mean, See Table IV 


t T Independent variable for frequency 
distributions, in standardized nota- 
tion, see equation (5-4) 


TANML Intermediate variable for computing 
YBAR of the Type Curve 


TEMPI A temporary storage cell for an 
intermediate value in the calculation 
of the Type IV curve 


TEDE Term for determining whether a data 
value is within overall interval limits 


o THETA . The angle equal to T(EU) + т/2 Гог 
use in the F(RV) calculation 


TOTAL Total number of property samples for 


a field 
TR Intermediate variable for computing 

YBAR of the Type Curve 

Ho U2 E moment about the mean, see 
equation (3-3) 

Hz 92 Bs moment about the mean, see 
equation (3-3) 

y U4 „EN moment about the mean, See 
equation (3-3) 

Ur ne D moment about the mean, see 
equation (3-2) 

¬ | ШЕЕ First moment of the histogram about 


the arbitrary chosen midpoint, See 
equation (3-6) 





Mathematical 
Notation 


1 
Ho 


Hz I 


73 


ғ. 


P4 


PA 


Fortran 


Notation 


UP(I) 


UMOD2 


UMODS 


UMODA 
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Description 


Second moment of the histogram about 
the arbitrary chosen midpoint 


Third moment of the histogram about 
the arbitrary chosen midpoint 


L'th moment of the data about the 
midpoint of the interval which con- 
tains the mean, see equation (3-1) 


ТЕ MOE O Nosvcnaerpedquenme 
distribution 


2'nd moment about the mean incorporat- 
ing Sheppard's corrections, see 
equation (3-10) 


3'!rd moment about the mean, incorporat- 
ing Sheppardlereoprreer sonen SEC 
equation (5-11) 


lith moment about the mean, incorporat- 
ing Sheppard's corrections, see 
equation (3-12) 


A parameter of the Pearson Type IV 
Curve, see equation (2-50 


A term in Bartlett's test equal to 
IMT 

i 
Degrees of freedom associated with each 
sample variance in Bartlett's test, 
see equation (5-13 


A parameter in the derivation of the 
constant for Type I Curve, see 
equation (2-26) 


The independent variable for the 
frequency distribution wherein each 
interval has unit width 


Equivalent of MP but in floating 
point notation 


A term used to compute 52 p’ see 
equation (5-6) 


The mean of a set of porosity data 





Mathematical Fortran 


Notation Notation 
Хы 
Xow 
XX 
XZERO 
Y Y 
YAVG 
YNORM 
2 Z 
2 
ZERO 
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Description 


The mean of a set of oil Saturation 
Gata 


The mean of a Set of water saturation 
data 


The values of the independent variable 
for a frequency distribution corresponds 
ing to the left, middle, and rignt 

ends of an interval. These intervals 
аге КГ туа 


The value of XX at the left extremity 
of a class interval 


The frequency for a Pearson type 
curve at any value of the independent 
variable, 5 


Integrated average value of frequency 
over a specific interval, see 
equation (3-15) 


Верона Ее ог frequency OM gic 
normal curve 


A parameter for the Pearson Type VI 
curve, see equation (2-36) 


The normal deviate corresponding to 
a desired confidence coefficient, 
see equation (6-6) 


A numerical constant, equal to zero 
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TECHNICAL AND STATISTICAL TERMS 


Term Notation Definition 
Array An arrangement of raw numerical data in 
ascending or descending order of magnitude 
Chi-Square хе А measure of the discrepancy existing 
between observed and expected frequencies 
Class Number of individuals belonging to each 
Frequency class in summarizing large masses of 


Degrees of 
Freedom 


Gamma 
lunes Lon 


Frequency 
Пра ет расло 


Histogram 


Mean 


Median 


Mode 


Moments 


Normal 
Distribution 


fe 


POS 


Ра! 


El 


raw data 


The number N of independent observations 
in the sample (i.e., the sample size 
minus The number K of population para- 
meters which must be estimated from 
sample observations 


A transcendental function used frequently 
in statistics and defined by equation NS 


A tabular arrangement of data by classes 
together with the corresponding class 
frequencies 


A set of rectangles having: (a) bases 
on a norizontal aris MC h) 
centers at the class marks and lengths 
equal to the class interval sizes, (b) 
areas proportional to class frequencies 


Arithmetic average of a set of data 


The middle value of a set of numbers 
arranged in order of magnitude 


That value of a set of numbers which 
occurs with the greatest frequency 


A term used in the measurement of 
dispersion of a distribution and defined 
by equation (5-2) 


Also referred to as Gaussian curve, 
normal curve of error, or normal 
probability distribution. A bell- 
shaped distribution defined by equa- 
tion (2-40) 





Term 
ou 
Saturation 


Permeability 


Population 


FOPOSTEN 


Range 


Skewness 


Standard 
Deviation 


Standardized 
Unite 


Standardized 
UNIE 


Unimodial 


Variance 


Water 
Saturation 


Notation 


Е 


O 


W 


Definition 151 


A measure of the oil present within a 
ре оце 
total fluid saturation within a rock 


A property of a porous medium and a 
measure of the capacity of the medium 
vo transmit fluids. Usually measured 
in dareiessor milliq o reres 


Term used in statistics to refer to the 
hypothetical complete enumeration of 
facts in-a particular field of study 


Ine ratio or thi Eyo dopa Cen nRa TOC 
to the buik volume of that rock 
multiplied by 100 to express it in 
percent 


The difference between the largest and 
smallest numbers in a set of numbers 


The degree of asymmetry, or departure 


irom symmetry; Or a distribution, see 


equation (2-13) 


The root mean square of the deviations 
from the mean 


A term to express distance from mean to 
midpoint of histogram intervals in units 
of standard deviations, see equation (3-4) 


A term to express the distance between 
the mean and another specified value of 
a frequency distribution as units of 
standard deviations 


For frequency curves, those curves 
that have only one mode 


The square of the standard deviation 
A measure of the water present within 


a rock expressed as a percentage of the 
total fluid saturation within a rock 
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APPENDIX B 


FORTRAN PROGRAMS AND FLOW CHARTS 
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650 FORTRAN PROGRAM TO 
CALCULATE A GOOD THEORETICAL 
FIT OF A PEARSON TYPE 
FREQUENCY CURVE FOR A GIVEN 
OBSERVED DISTRIBUTION. 
THE PROCEDURE TO BE FOLLOWED 
1 TABULATE THE DATA USING 
CONVENIENT CLASS INTERVALS, 
2 CALCULATE THE MOMENTS ABOUT 
A SELFCTED VERTICAL». 
3 TRANSFER THE MOMENTS TO 
THE MEAN. 
4 APPLY SHEPPARDS CORRECTIONS 
TO THE MOMENTS, 
5 CALCULATE ALP3S ALP4 DELTA 
6 LOCATE THE MEAN (DAVG) 
7 DETERMINE FROM THE CHART 
WHAT TYPE OF CURVE TO USE 
PROGRAM BY J S VAN SCOYOC 
UNDER THE DIRECTION OF DR. 
FLOYD PRESTON 


DIMENSION DD( 100) +.G(100) »UP(8) 
sSUMXG(8)9X(100) »XG(8) 

READ1» INDEX »NDEX 

READ1» DMIN» DMAXs DELD»sDECM 


DETERMINE NUMBER OF INTERVALS 


AN =(DMAX-DMIN) /(DELD+DECM) 
IF (XCONF(1)) 42250950 
РОМСН1» АМ 

L = XFIXF (AN) 


DETERMINE CENTER OF INTERVALS 


DD(1) = DMIN +DELD/2.0 

IF (ХСОМЕ (1) ) 62 +70+ 70 
РОМСН1• DD(1) 

0 80 I = ls L 

DD( 1+ 1) =рО(1) + ОЕ О+1» 
ТЕ(ХСОМЕ(1)) 82,90,90 
РОМСН1» 00 


COUNT DsS IN EACH INTERVAL 
AND COUNT TOTAL D»S AND SUM Ds 


AK=008 

SUMD=0-6 

II = 1 

J=l 

А./-0. 

GO TO (150,180), INDEX 


READ IN DATA 
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150 
READ 1, IDs 


I PERM 


ISAMPs IDEPTs 


CONVERT TO FLOATING POINT 


DATA = FLOTF (IPERM) 

GO TO 270 

READ IN DATA 

READ1s IDs ISAMPs IDEPTs IPHIs = 
ISOsISW 


GO TO( 20092203240) sNDEX 


CONVERT TO FLOATING POINT 


DATA = ГЕОТЕ "ГЕТ 
COMO) 70 

DATA = FLOTF (ISO) 
GO TO 270 

DATA = FLOTF ( ISW) 


IF (DATA=DMIN) 140 5280 5280 

GO TO(2909330)911 

II = 2 

TEST = DMIN + DELD+DECM 

IF (DATA-TEST) 324093809380 
SUMD = SUMD + DATA 

Ad = Al Ш 

ЈЕ (ХСОМЕ (1) ) 362 + 37 0+ 370 
PUNCHI$SUMDsSDATAsSTESTsSA JS AKs Jes - 
І 

GO TO 140 

AK = AK + AJ 

G(J) = AJ 

AJ=0 

а-ы A 

ТЕЗ = ТЕ + DELD+L+ 

IF (TEST -IDMAX+DECM))330s»330s 
430 


CALCULATE MEAN 


DAVG = SUMD/ AK 
ТЕ ( ХСОМЕ (1})432.440.440 
PUNCH1IsSUMD sAK sDAVGsG 


DETERMINE INTERVAL WHICH 
CONTAINS MEAN 

DMP = (DAVG—-DMIN) 
MP = XFIXF (DMP) 
XMP = FLOTF (MP) 
IF (XCONF (1) 146224803480 
РОМСН1» ОМР»МР»ХМР 


/ (DEL D+DECM + 
sal 


PLACE INTERVALS ON UNIT BASIS 


DO 490 I = 19» 8 
SUMXG(I) = O, 
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NA=L+1 
DO 560 1 = 1+МА 
Xu) =WFLOTF TI)7-XMP 


IF (XCONF (1))512:520520 
PUNCH] o X(T) 


CALCULATE XG TO X8G 


ala meno Ae 

DO 540 J z ls. 

XG(J+1) = XG(J) * X(I) 
ІЕ(ХСОМҒ(1)1542»550»550 
РОМСН1» XG 

DO 560 J = 1, 8 


SUM XG TO X8G 


SUMXG(J) = SUMWKG(J) + XG( J) + 
IF(XCONF(1)) 56255705570 
PUNCH1s SUMXG 


CALCULATE UPRIMEIS 


DO 580 I = ls 8 

UP(I) = SUMXG (I) / AK 
IF (XCONF(1)) 58295909590 
PUNCH1s UP 


CALCULATE U2 TO U8 


U2=UP(2)-UP(1)*UP(1) 
U3=UP(3)-3.*UP(1)*UP(2)+20+* - 
UP(1)*UP(1)*UP(1) 

UAzZzUP(4)-A4. *UP(1)*UP( 3) 46. X - 
UP(CI)*UP(1)*UP(2)-3. *UP( 1)0* - 
UP(1)*UP(1)*UP(1) 

SIGMA=U?#*,5 

IF(XCONF(1)) 6225623056230 
PUNCH1sU1sU2sU3,U4 


APPLY SHEPPARDS CORRECTIONS 


UMOD2=U2—-le/1l2e 

UMOD3=U3 
UMOD4=U4-0.5*U2+7./2408 
IF(XCONF(1)) 65256605660 
PUNCH1 »s UMOD2 s UMOD3 »UMOD4 


CALCULATE VALUES NEEDED 


ALP 3S=UMOD3*UMOD3/ (UMOD2*UMOD2 = 
*UMOD2) 

ALP3=ALP3S##,5 

ALP4=UMOD4/ (UMOD2#UMOD2) 
DELTA=(2.#ALP4-3.#ALP3S-6.)/ + 
(ALP4+3.) 
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00790 I=1»NA 
РОИМСН] » То00 (1) »Х (1) +6(1) 
CONTINUE 
PUNCH1 sALP3 sALP3SsDELTAsSIGMAs 
ОРУТ 
GO TO 30 

END 
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156 
FLOW CHART TO IDENTIFY THE PEARSON TYPE CURVE 


CORRESPONDING TO A GIVEN DATA SET 


READ, | 
INDEX, NDEX, DMIN, DMAX, DELD, DECM 
DD(1) = DMIN + DELD/2 


І = І 


DD(I-1) - DD(I) 4 DELD + 1 


AK 
АЈ 


Тр, ТЗАМР, ТРЕРТ |. 1 2 ID, ISAMP, IDEPT 
IPERM IPHI, ISO, ISW 


He 


DATA = ISO 


C4 Hl 
| | 
HH 







II 22 TEST = DMIN + DELD + DECM 


6% 
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AK = AK +d 


il 
I~ 


SUMD = SUMD + DATA DATA: TEST 


б DOES + 1 


(С = TEST: DMAX + DECM 
- | 


= Ј + 1 
DAVG, MP, XMP, DMP 
Tel) TEST = TEST + DELD + 1 
gc 


SUMXG (I) = | -——— ИЕ 


2 
xo (141) exe(2)* X (1) 
UP(I)=SUIEKG (I) /AK 
2 
SIGMA,U2,U3,U4 UMOD2, UMOD3 , UMOD4 ALP3, ALP4, DELTA 
PUNCH 
NA, N, DMIN, DMAX, DELD, DAVG NR Ае 
Ga renes 7, рыт), х(т), а(т) 


PUNCH 2 zc NA 
ALP3, ALP35, DELTA, SIGMA, UP(1 


8) 
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PEARS@GM FIT РОК TYPES 193% 496 
DATA FOR READ2,READ3, AND 
READ4s ARE SUPPLIED BY THE ` 
IDENTIFICATION 

NUMBER (FOR READ5»s )=NUMBER 


OF DIFFERENT TYPES OF CURVES + 
WHICH ARE TO BE FITTED TO THE 
SAME DATA 


NTYPE=CURVE TYPE NUMBER 

NTYPE=1 FOR PEARSON TYPE 1 
NTYPE=2 FOR PEARSON TYPE 3 
NTYPE=3 FOR PEARSON TYPE 4 
NTYPE=4 FOR PEARSON TYPE 6 


THIS PROGRAM NEEDS SINF AND * 
XCONF SUBROUTINES AND SPECIAL + 
GAMMA FUNCTION ROUTINE 


DIMENSION DD(50)+X(50)9G[(50)» 
ХХ (3) э\ (3) УМОВМ (3) › СОМ (5) » 
CST(16),EU(16),EZ(16)» 

ТЕЕТА 167 

ZERO= 0% 


READ1 sCSTsEU 

DO 501 I=1+16 
THETA(1)=3.1415927*EU(1)+ 
15707963 
EZ(1)=SINF(THETA(I)) 


READ2 sNAsTOTAL sDMINsDMAXsDELD > + 
DAVG 

AN=FLOTF (NA) 

DOAI=1»NA 


READ3esJeDD( 1) o»X(I) sG(1) 
READ4 + ALP29ALP3S+DELTA9+SIGMA - 


ЦР] 

О= АЫ РЗжА РЗ - 4 • ХрЕ| ТА # 
(DELTA+20) 

1F(D)89999 

D=-D 

DSQRT=D**.5 


READS sNUMBR 
J=] 


READ6»NTYPE 
00141=1»5 
CON(1)=08 
PUNCH1»NTYPE 
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Om O O О 


Q O (O O O O O O O Or O O O O OO O O OLO = O QOO О 


ON OO 


ооомооооооо-оо-о 


NPCHz417 

ЈУРИ Хе (1 ))417• 17917 

РОМСНТ] о МРСНОМТУРЕ МОМВЕ ФМА » 0 o 
DSQRT AN 

GO T0( 1833059365318) »sNTYPE 


EM1=—((1e+DELTA)/DELTA)*¥(le- 
Да ДО ОК Те 

СОМ (1 } =ЕМ1 
ЕМ2=—((1•+0рЕ ТА) /РЕ ТА) #( 1 + 
ALP3/DSQRT)-le 

CON(2)=EM2 
R1=(-ALP3+DSORT)/(2e*DELTA) 
CON(3)=R1 
R2=(-ALP3-DSORT)/(2¢*DELTA) 
CON(4)=R2 


GO T0(27929929548) sNTYPE 
С=ТОТАС/ ( $1СМАХ (К2-К1)**(ЕМ] + 
ЕМ2 +1.) %*САМАЕ (ЕМ1 +] ) *GAMAF (EM 
2+1.) ) 

СОМ (5) =С 

GO. .T0. 541 

STOP 


A=2e/ALP3 

CON(1)=A 

ASQ=A#A 
CsA**(ASQ)/USIGMA*EXPEF(ASQ) * 
GAMAF(ASQ) ) 

CON(2)=C 

Comoe) 


R=ALP3/(2.*DELTA) 

CON(1)=R 

EM=1+./DELTA+20 

CON(2)=EM 

S=DSQRT/(2e¢*DELTA) 

CON(3)=S 

P1=2.*EM-20 

V=-2.* (1. +DELTA)*XALP3/(DELTA* 


DSORT) 
CON(4)=V 
ТЕМР 1 =2. *ЕМ-1 
50М=0, 


рО 502 1=1,16 

Vata (G4 VER) EXPER IVETHETACT) ) 
SUM=SUM+Y*CST (1) 
FRV=EXPEF(-1.5707963*V ) *SUM+ 
301415927 
C=TOTAL*S#**TEMP1/(SIGMA#FRV) 
CON(5)=C 

COMO 51 


ALPHA=R1-R2 

C=TOTAL*GAMAF (-EM2)/(SIGMAX* 
GAMAF(EM1+1.)%*GAMAF(-EM2-EM1- 
1.) *XALPHA**(EM1+EM2+10)) 
CON(5)=C 
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CHISQzO, 

NPCHz452 

Tr CACORE U1 52» 5255 2 
РОМСН1»МРСН»СОМ 

BEL SQzO, 
XZEROzX(1)-.45 
PROP=DMIN+DELD*#.5 


T=] 
XX(1)=XZERO 
ХХ(2)=Х (1) 

ХХ (3) ЕХ2ЕКО+] в 
L=1 

XBAR=XX(L) 


T= (XBAR-UP1)/SIGMA 

УМОКМ (|) =. 398942] *ЕХРЕР (-Т*®Т/2 = 
e)* TOTAL/SI GMA 

GO ТО(63 968.70 ,65 ) »sNTYPE 


YBARZ(C*(T-RI)**EMI)*((R2-T)** 
FM?) 
So RER] 


Z=(T=R2) 
YBARzC*Z**FEM2*(Z-ALPHA)?**EMI = 
GOTO 91 


YBARZCX*(A-T)**(ASQ-14 ) *EXPEF * 
(-А%Т) 
GOTO 


TR=(T+R) 

ARG=TR/S 

XSQ=ARG*ARG 

COE=3. 
IF(XSQ-1.)83»175»76 
ТАММ] =. 78539816 

GO TO 89 
TANM1=1.5707963-1./ARG 
Qzle 

00811=1%+9 
TANMIZTANMI-4Q/(ARG*XSQ*COE) * 
XSQzXSQ*ARG*ARG 
COE=COE+2¢ 

Q=-Q 

Go TO 89S 

TANM1=ARG 

Q=-le 

DO 881=199 
TANM1=TANM1+Q*ARG#XSQ/CCE 
XSQzXSQ*XARG*ARG 
COEZCOE-*2. 

Ос-0 

АВС1 =\* (1.5707963-ТАММ1 ) 
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PUNCHI F NPCHs I IsL + XX + YBAR +T, 
ҮМОКМ»Х5О»ТАММ1 
IF(L-3)93,95,95 

L-L-*1 

GO TO 59 


ҮМАМС= ( ҮМОКМ ( 1 ) +ҮМОКМ (3 ) +4, % + 
YNORM(2))/30 

YNORM(2))/6. 

YAVG=[(Y(1)+Y(3)+4.*Y(2))/60 + 
РОМСН191»Х( І) «РКОР»С( І)»Ү(2)» 
YAVG+YNORM(2) 

IF(I-NA)99,1115111 


BSQ=G(1)-YNAVG 
BSQ=BSQ*BSQ 

IF(YNAVG) 10151025101 
BELSQ=RBELSQ+BSQ/YNAVG 
CHI =G(1)-YAVG 

CHT ЕНИ Xx CET 

IF(YAVG) 10441055104 
CHISQ=CHISQ+HCHI/YAVG 
XZERO=XZERO+1 « 

Y(1)=Y(3) 
YNORM(1)=YNORM(3) 
PROP=PROP+DELD 

МРСН=4109 

ТЕ (ХСОМЕ (1) ) 4109» 109» 109 
РОМСН1»МРСН» І91»СНІ»СНІ50»В50» 4 
BELSQ 

1=1+1 

G0 ©0 55 


PUNCH1sZERO 
PUNCH1sDMINsDMAX sDELD sDAVGs 
SIGMAsALP3sDELTA 
РОМСН1»2ЕКО 
РОМСН1» СОМ» СНІ50»ВЕІ 50 

IF (J-NUMBR) 11651219121 
J=J+1 

PAUSE 


GO TOLLIS 
PAUSE 

GU TOTO 
END 
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16 
GENERALIZED FLOW CHART FOR CURVE FIT TO E 


PEARSON TYPE DISTRIBUTION 


READ 
TABLE FOR F(RV) 


NA, TOTAL, DMIN, DMAX, DELD, DAVG 
S DDU sS 


READ 
ALP3, ALP3S, DELTA, SIGME Al 


READ 
NIYPE 


CALCULATE NECESSARY FIXED 
PARAMETERS AND CONSTANTS C 
FOR SELECTED NTYPE 


CALCULATE PROP, Y(2), 
YAVG, YNORM(2), FOR 
.EACH INTERVAL I 








O | 105 


PUNCH 
I, X(I), PROP, G(I), Y(2), YAVG, YNORM(2) 


2 
ІМ 


CALCULATE CHI AND BSQ 
FOR EACH INTERVAL AND 

SUM TO GIVE CHISQ 
AND BELSQ 






= 


DMIN, DMAX, DELD, DAVG, 


9 SIGMA, ALP3, DELTA 


PUNCH 
CON, CHISQ, BELSQ 





FIGURE 32 
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EXPLANATION OF FLOW CHART SYMBOLS 


- > || 


A Read or Punch instruction. 
Data are read or punched in the 
согабы убал lel CO EIEN 


An arithmetic calculation. 


A test of the difference A - B 

Tor negativity; рост олур 
zero. The program branches to 

the next appropriate Calci ае ой 
sequence depending upon the 
numerical value of this difference. 


A multi-way branch depending 
upon the numerical value of 
the inclosed index. 


A jump is made in the flow 
chart to the circle with the 
Same enclosed letter. 


FIGURE 55 
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FIGURE 34 MAP OF WELLS CORED--FIELD 1 
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FIGURE 35 MAP OF WELLS CORED--FIELD 2 
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FIGURE 36 MAP OF WELLS CORED--FIELD 3 
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